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The Hubble diagram (HD) is a plot of the measured distance modulus versus 
the measured redshift, with the slope giving the expansion history of our uni- 
verse. In the late 1990's, observations of supernova out to a redshift of near unity 
demonstrated that the universal expansion is now accelerating, and this was the 
first real evidence for the mysterious energy we now call Dark Energy. One of 
the few ways to measure the properties of Dark Energy is to extend the HD to 
\ higher redshifts. Many models have been proposed that make specific predictions 

as to the shape of the HD and so this offers a means of testing and eliminating 
models. Taking the HD to high redshifts provides a way to test models where 
their differences are large. For example, in a comparison of the now concordance 
model (a flat universe with Q>m — 0-27 or so and constant Cosmological Con- 
stant) with a representative model of evolving Dark Energy (for which I'll take 
w(z) = —1.31 + 1.482 from Riess' analysis of the gold sample of supernovae), the 
predicted distance moduli differ by 0.15 mag at z = 1.7 and 1.00 mag at z — 6.6. 
The only way to extend the HD to high redshift is to use Gamma-Ray Bursts 
(GRBs). GRBs have been found to be reasonably good standard candles in the 
usual sense that light curve and/or spectral properties are correlated to the lu- 
minosity, exactly as for Cepheids and supernovae, then simple measurements can 
be used to infer their luminosities and hence distances. GRBs have at least five 
properties (their spectral lag, variability, spectral peak photon energy, time of the 
jet break, and the minimum rise time) which have correlations to the luminosity 
of varying quality. All of these properties provide independent distance informa- 
tion and their derived distances should be combined as a weighted average to get 
the best value. For GRBs which have an independently measured redshift from 
optical spectroscopy, we have enough information to plot the burst onto the HD. 
In this paper, I construct a GRB HD with 69 GRBs over a redshift range of 0.17 
to >6, with half the bursts having a redshift larger than 1.7. This paper uses 
over 3.6 times as many GRBs and 12.7 times as many luminosity indicators as 
any previous GRB HD work. For constructing the GRB HD, it is important to 
perform the calibration of the luminosity relations for every separate cosmology 
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considered, so that we are really performing a simultaneous fit to the luminos- 
ity relations plus the cosmological model. I have made detailed calculations of 
the gravitational lensing and Malmquist biases, including the effect of lensing 
de/magnification, volume effects, evolution of GRB number densities, the GRB 
luminosity function, and the discovery efficiencies as a function of brightness. 
From this, I find that the biases are small, with an average of 0.03 mag and 
an RMS scatter of 0.14 mag in the distance modulus. This surprising situation 
arises from two causes, the first being that burst peak fluxes above threshold 
do not vary with redshift and the second being that the four competing effects 
nearly cancel out for most GRBs. The GRB HD is well-behaved and nicely delin- 
eates the shape of the HD. The reduced chi-square for the fit to the concordance 
model is 1.05 and the RMS scatter about the concordance model is 0.65 mag. 
This accuracy is just a factor of 2.0 times that gotten for the same measure from 
all the big supernova surveys. I claim that GRBs will not suffer from effects due 
to evolution in the progenitors as we look back in time, with the reason being 
that the luminosity indicators are results of light travel time delays, conserva- 
tion of energy in the shocked material, and the degree of relativistic beaming, 
with these not changing with metallicity or age. That is, even though distant 
bursts might be more luminous on average than nearby bursts, the luminosity 
indicators will still operate to return the correct luminosity. I fit the GRB HD 
to a variety of models, including where the Dark Energy has its equation of state 
parameter varying as w(z) = Wo + w a z/(l + z). I find that the concordance model 
is consistent with the data. That is, the Dark Energy can be described well as a 
Cosmological Constant that does not change with time. 

Subject headings: Gamma-Ray: Bursts - Cosmology: cosmological parameters - 
cosmology: observations - cosmology: dark matter 

1. Introduction 

The Hubble diagram (HD) is a plot of distance versus redshift, with the slope giving the 
expansion history of our Universe. This expansion history depends on the amount of mass 
in our Universe (both normal and dark matter) as well as the Dark Energy. So by making 
observations of distances over a wide range of redshifts, we can measure the expansion history 
and place significant constraints on models of the Universe. Indeed, it was this program that 
lead to the discovery of Dark Energy (Perlmutter et al. 1999), while observational knowledge 
about Dark Energy arises solely from the constraints on the expansion history. The program 
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of measuring the HD with supernovae has now occupied a vast amount of telescope time 
during the last decade and it will likely lead to a dedicated satellite over the next decade. 

Supernovae are now seen as near-ideal 'standard candles' for purposes of the Hubble di- 
agram. (The term 'standard candle' is commonly used to indicate a light source by which the 
luminosity can be determined regardless of whether the objects in that class all have the same 
luminosity. For example, both Cepheids and Type la supernovae are called standard candles 
even though they vary in absolute magnitude by up to four magnitudes.) But it is worthwhile 
remembering that only a decade ago the situation was that only one moderate-redshift event 
had ever been seen (N0rgaard-Nielsen et al. 1989) and there was still widespread debate as to 
whether supernovae even were 'standard candles' (van den Bergh 1993; 1996; van den Bergh 
& Pazder 1992; Phillips 1993; Tamman & Sandage 1995). This situation changed only when 
Hamuy et al. (1996) used a large and well-observed data set to prove that Type la supernovae 
were indeed 'standard candles', Perlmutter et al. (1995) found efficient methods to discover 
distant events, Kim, Goobar, & Perlmutter (1996) solved the k-correction problem, and Perl- 
mutter et al. (1997) started reporting on high-redshift light curves. As is often the case with 
large and complex programs, the preliminary results (Perlmutter et al. 1997) gave a com- 
pletely different conclusion than the more settled later conclusions (Perlmutter et al. 1998; 
1999). The nearly simultaneous results from an independent group of supernova hunters 
(Riess et al. 1998) provided the community with a good sense about the validity of the 
conclusions. Later massive campaigns involving ground-based (Astier et al. 2006) and HST 
observations (Riess et al. 2004) have confirmed the earlier results. Nevertheless, as is usual 
for important claims, a variety of troubles were raised, including grey-dust (Aguirre 1999a; 
1999b), evolution of the supernova progenitors (Dommguez, Hoflich, & Straniero 2001), and 
refraction by Lyman-a clouds (Schild & Dekker 2005). But what really convinced the as- 
tronomy and physics community is that other independent lines of evidence have now been 
made that confirm the cosmological parameters derived from supernovae (Lineweaver 1999; 
Spergel et al. 2003; Eisenstein et al. 2006). These results are consistent with a cosmology 
where the Universe is flat with VL M = 0.27 and Ha = 0.73, in which the Universe has an 
age of close to 14 billion years and will expand forever. The observed acceleration of the 
universal expansion is attributed to a mysterious energy labeled 'Dark Energy'. 

The best way to measure properties of the Dark Energy seems to be to measure the 
expansion history of our Universe. To this end, the SNAP satellite (Scholl et al. 2004) has 
been proposed to determine the distances and redshifts of two thousand supernovae per year 
out to redshift 1.7 with exquisite accuracy. The default expectation is the simplest model 
for the Dark Energy, where it does not change in time. This can be parameterized with 
the equation of state of the Dark Energy behaving as p = wpc 2 , where p is the pressure, p 
is the density, c is the speed of light, and w is a dimensionless constant that might change 
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with time. The concordance case has w — — 1 at all times, and this is the expectation of 
Einstein's cosmological constant, or if the Dark Energy arises from vacuum energy. Given 
the strong results from supernovae for redshifts of less than 1, the frontier has now been 
pushed to asking the question of whether the value of w changes with time (and redshift). 
Lacking any physical understanding of such changes, we can parameterize the changes with 
low-order expansions such as w(z) = wo + w'z or w(z) = wo + w a z(l + z)^ 1 (Linder 2003). 
So far the Supernova Legacy Survey results (Astier et al. 2006) have not been used to check 
the constancy of w, but the 'gold sample' of supernovae (Riess et al. 2004) shows that the 
best fit has w = —1.31 and w' = 1.48. 

If the Dark Energy changes with redshift, the best plan is to measure it over a wide 
range of redshifts, but supernovae cannot be detected out past redshift 1.7 even with the 
next generation of satellites. Many alternatives to the concordance (w = —1) case have 
been proposed, including Weyl gravity (Mannheim 2006) and many versions of quintessence 
(Szydlowski, Kurek, & Krawiec 2006), and at this point there is no strong reason to prefer 
the concordance. The concordance cosmology makes a unique prediction as to the shape 
of the HD out to high redshift whereas alternatives make distinct and different predictions, 
so here is a perfect case where observations can test models of the Dark Energy. After 
all, the concordance cosmology has only been around for a few years, whereas the previous 
concordance case (with Q M = 1 and JIa = 0) was held just as dearly until its demise. So a 
strong imperative in the quest for Dark Energy is to extend the HD to high redshift. 

Gamma-Ray Bursts (GRBs) offer a means to extend the Hubble diagram to at least 
redshifts of >6. The reason is that GRBs are visible across much larger distances than 
supernovae. For example, the average redshift of GRBs detected by the Swift spacecraft 
(Gehrels et al. 2004) is 2.8, while the range is from near zero out to >6 (Jakobsson et al. 
2006). GRBs are now known to have several light curve and spectral properties from which 
the luminosity of the burst can be calculated (once calibrated), and these make GRBs into 
'standard candles'. Just as with supernovae, the idea is to measure the luminosity indicators, 
deduce the source luminosity, measure the observed flux, and then use the inverse-square 
law to derive the luminosity distance for plotting on a HD. To be placed on the HD, each 
GRB must have an independently measured redshift, usually from absorption lines in the 
optical spectra of the GRB afterglow or from optical spectra of the host galaxy. To date, 
roughly 95 GRBs have a measured redshift (Greiner 2006). 

I presented the first GRB Hubble diagram in three public talks in March and April 2001 
based on 8 GRBs (expanded to 9 GRBs in Schaefer 2001). The first published GRB HD 
appeared in Schaefer (2003a) based on 9 GRBs using two luminosity indicators. Half a year 
later, another GRB HD appeared with 16 bursts based on the independent result that the 
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gamma-ray energy of the bursts is a constant after correction for the beaming angle (Bloom 
et al. 2003). Since then, various groups have used one single luminosity indicator each to 
construct HDs or to similarly constrain cosmology (see Table 1). This table shows that all 
prior work has been made with only a small fraction of the available data; both taking only a 
fraction of the useable bursts and in using only one or two luminosity indicators while more 
information from other indicators is ignored. All the prior GRB HDs have had too large 
error bars to provide useful constraints on cosmology. 

The purpose of this paper is to make use of all possible GRBs and to make simultaneous 
use of all possible luminosity indicators. As such, this paper will use over one order of 
magnitude more luminosity indicators than any previous work (see Table 1). A preliminary 
version of this work was presented at the January 2006 meeting of the American Astronomical 
Society (Schaefer 2005), where I reported that the first results based on 52 GRBs showed 
a 2.5-sigma deviation from the unique prediction of the concordance cosmology. Now, with 
17 new GRBs plus substantial amounts of new data for the original 52 events as well as 
with some small improvements in analysis method (see below), the conclusions presented 
in this paper have a somewhat smaller deviation than in my preliminary results. Just as 
with the change between the preliminary and final results of Perlmutter et al. (1997; 1999), 
I expect that the results in this paper will be closely checked by independent analysts, 
duplicated with independent events (collected over the next two years by Swift), and tested 
for potential problems. A lesson learned from the supernova cosmology experience is that 
any results from the HD for 1.7 < z < 6.6 will only be accepted by the community after 
independent methods reproduce the same conclusions. 



2. Observations 

For a GRB to be placed on the HD, it must have a measured brightness, measured 
luminosity indicators, and an independent redshift. The brightness and luminosity indicators 
require that a light curve and/or a spectral fit be available, while the redshift requires an 
optical spectrum of the afterglow or the host galaxy. Greiner (2006) keeps an updated web 
page with all the bursts and which ones have reported redshifts, and this is the primary 
basis for selecting bursts. This section will tell about the bursts selected and their measured 
properties. 

Out of the roughly 95 GRBs with reported redshifts, I have identified 69 to-date (with 
a cutoff at 7 June 2006) as having adequate data for placement on the GRB HD. These are 
listed in Table 2. 
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Some bursts with reported redshifts are not included for a variety of reasons: GRBs 
050509, 050709, 050724, 051221, and 060502B are all short-hard bursts that belong to a 
different population than normal GRBs, and the luminosity indicators do not apply to this 
separate class of bursts. GRB 051227 might be a short-hard burst and it only has a limit 
on its redshift. GRBs 980329 and 011030 only have limits on their redshifts. GRBs 980326, 
011121, 020305, 050803, and 060123 have redshifts of too low a confidence to be used. GRBs 
050730, 050814, 051016, 060218, 060512, and 060522 are all too faint to have any useful 
constraints on their luminosity indicators. GRBs 000214, 000418, 000301C, and 031203 do 
not have adequate light curve or spectral data available. GRBs 980425, 020819, and 050315 
are outliers of various types and will be handled separately in section 3. 

For each selected GRB, I have tabulated a range of primary information in Table 2. 
Column 1 gives the six-digit GRB identification number in the usual format of YYMMDD. 
Column 2 identifies the satellite experiment that is the source of the light curve used in each 
case. Column 3 gives the redshift (z) of the GRB rounded to the nearest 0.01. Column 4 gives 
the reference number (see the footnote to the table) for the source of the redshift. Column 5 
gives the observed peak flux (P) and its one-sigma uncertainty over the brightest one-second 
time interval during the burst. Some uncertainties are not reported in the literature, so I have 
adopted a value of P/10 as the error bars for P as indicated by values in square brackets. The 
units are either photons/cm 2 /s for values larger than 0.01 or in units of erg/cm 2 /s for values 
smaller than 0.01 (as expressed in scientific notation). All but four of the GRBs have the 
quoted P values for a time interval length of 1 second (or 1.024 seconds), and the exceptions 
are marked with a footnote that also gives a conversion factor to convert to a 1 second time 
interval as based on direct calculation from the observed light curves. Column 6 gives the 
quoted energy range for the given P values, E m i n and E max . Column 7 gives the reference 
for the quoted peak flux data. Column 8 gives the peak energy of the vF v spectrum (E peak ). 
The error bars are often asymmetric and so I first give the uncertainty in the direction of 
lower E peak followed by the uncertainty in the direction of higher E peak . Again, one-sigma 
error bars are often not quoted in the literature and so I have adopted the typical values as 
shown in square brackets. For some GRBs, the E peak is not reported and so I have adopted 
an approximate E peak value based on the best fit luminosity relation for E peak along with 
very large adopted error bars. It is important to note that these adopted E peak values are 
only used for extrapolating the spectra to the bolometric value and this assumption has 
little effect because the observations cover the energies with most of the flux. Columns 9 
and 10 give the values for the low-energy and high-energy power law spectral indices (a and 
j3 respectively) along with their one-sigma uncertainties. When a value of a is not known I 
adopt a value of —1.1 ± 0.4. When a value of (3 is not known I adopt a value of —2.2 ± 0.4. 
These spectral parameters are only used for estimating the bolometric peak fluxes and the 
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adopted values make relatively little difference while the uncertainties are propagated into 
the bolo metric peak fluxes. Column 11 gives the reference for the spectral information. 

Some GRBs have a measured jet break time (tj et ) when the afterglow brightness has 
a power law decline that suddenly steepens due to the slowing down of the jet until the 
relativistic beaming roughly equals the jet opening angle, 9j et (Rhoads 1997). The measured 
tj et value can be used to deduce 9j et and hence convert the 'isotropic' energy (E ltiso ) into the 
total energy in gamma-rays emitted by the burst (E y ). From Sari, Piran, & Halpern (1999), 

9 jet = 0.161[t jet /(l + z)f 8 (n v,/E m52 ) 1/8 , (1) 

where z is the redshift, tj e t is the jet break time measured in days, n is the density of the 
circumburst medium in particles per cubic centimeter, ?7 7 is the radiative efficiency, and 
E lt iso,52 is the isotropic energy in units of 10 52 erg for only an Earth-facing jet. Here, 9j et 
is in degrees and it is the angular radius (i.e., the half opening angle) subtended by the jet. 
The isotropic energy is calculated as 

E^ iso = 47cd 2 L S bolo {l + z)- 1 , (2) 

where is the luminosity distance of the burst and S bo i is the bolometric fluence of the 
gamma-rays in the burst. The total collimation-corrected energy of the GRB is then 

E 1 = (1 - cosOjet) E J>iso , (3) 

where the beaming factor, F beam is 1 — cos9j et . 

One of the luminosity relations connects E 1 with E pea k, so bursts with an observed jet 
break can be used as a luminosity indicator. Unfortunately, only a little less than half of the 
bursts have observed jet breaks. The important observational quantities have been listed 
in Table 3 for these bursts. Column 1 gives the GRB name. Column 2 lists the observed 
fluence (S) and its uncertainty. Some error bars are not reported in the literature, so I have 
adopted a conservative 10% error as indicated by square brackets. All values are in units of 
erg/cm 2 . Column 3 gives the energy range (in units of keV) for the reported fluence, E min 
and E max . Column 4 gives a reference for the source of the quoted fluence. Column 5 gives 
the bolometric fluence and its uncertainty (in units of erg cm -2 ) as discussed in the next 
paragraph. Column 6 lists the measured time of the jet break (tj e t) in units of days, while 
column 7 gives the reference for these values. Columns 8 and 9 give the derived 8j e t values 
(in degrees) and the reference. For many bursts, the value of 9j et is based on a detailed 
model fitting of light curves in many energy bands (e.g., Panaitescu & Kumar 2002). In the 
absence of these detailed fits, it is reasonable to adopt ?? 7 = 0.2 and n = 3 cm -3 as these 
values are present in the equation for 9j et to the 1/8 power. Finally, column 10 gives the 
derived value for F beam and its uncertainty. 
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The reported brightnesses (peak fluxes and fluences) are given over a wide variety of 
observed band passes, and with the wide range of redshifts these correspond to an even wider 
range of energy bands in the rest frame of the GRB. One way to put all these brightnesses 
onto a consistent basis is to derive the bolometric brightness, where the measured spectrum 
is extrapolated to high and low energies and integrated over all energies. This can be done 
with fair accuracy for GRBs as the bulk of the energy comes out in and near the observed 
band passes so that uncertainties in the extrapolation are generally small. In practice, I will 
integrate under the measured spectral shape from photon energies of 1 to 10,000 keV in the 
rest frame of the GRB. I adopt the shape of the GRB spectrum to be a smoothly broken 
power law (Band et al. 1993) as 

/ A E* e-( 2 +") £ /*W if E < [(a - 0)/(2 + a)]E peak 
HE) = \BEP otherwise (4) 

Here $ is the usual differential photon spectrum (dN/ dE) as a function of the photon energy 
(E), Ep ea k is the photon energy at which the vF v spectrum is brightest, a is the asymptotic 
power law index for photon energies below the break, and f3 is the power law index for photon 
energies above the break. The two normalization parameters (A and B) are chosen to ensure 
continuity at the break and to match the observed burst brightness. If the observed peak 
flux (P) is reported in units of erg/cm 2 /s, then the bolometric peak flux will be 



/•10 4 /(l+z) pE max 

Pbolo = P / E$dE I \ E§dE (5) 

A/(l+z) JE mi „ 



If the units of P are photon/cm 2 /s, then 

r-10 4 /(l + 2) pE. 



Pbolo — P 

'1/(1+*) 



plU-/(L-fZ) rJ-^max 

/ E<$>dE I I <$>dE (6) 



The result will be the bolometric peak flux over a one second interval with units of erg/cm 2 /s. 
Similarly, for the observed fluences (S all in units of erg/cm 2 ), we can calculate the bolometric 
fluence as 

rl0 4 /(l+z) rE max 

S bolo = S / E$dE I \ EQdE. (7) 

The uncertainties on both Pbolo and Sbolo are both calculated by propagating the uncertainties 
in P, S, a, and (5. The resulting bolometric peak fluxes and fluences are tabulated in Table 
4. The bolometric fluences are given in Table 4 only for those bursts with known jet break 
times, as that is all that is needed for the purposes of this paper. The isotropic luminosity 
is given by 

L = Arrd 2 L P b oi . (8) 
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To calculate either L or E~ t (for use in the luminosity relations), we have to know the 
luminosity distance (c?l), which depends on the cosmological model and the redshift. 

The lag (tj 09 ) of a GRB is the time shift between the hard and soft light curves, with the 
soft photons coming somewhat later than the hard photons. Schematically, we can view this 
as the delay in the time of peaks in the, say, 100-300 keV versus 25-50 keV energy bands. But 
in practice, only the brightest bursts have their times of peaks defined well enough to make 
this a useful definition. For most bursts, the normal Poisson noise in the light curve will 
create substantial errors if only the bins with the highest fluxes are compared. A practical 
definition of the lag can be taken as the offset that produces the maximum cross correlation 
between the hard and soft light curves. And for faint bursts, the cross correlation itself will 
have substantial noise so that it is best to fit a parabola around the maximum. To avoid 
the noise added by time intervals when the light curve is dominated by ordinary background 
noise, the cross correlation is performed only for times when the total light curve (summed 
over all available energy bands) is brighter than 10% of the peak flux in that light curve. The 
choice of energy bands for the hard and soft light curves should be such that there is as wide 
a separation as possible yet with the high energy band still having significant flux. A good 
choice is the original choice for the BATSE data (Band 1997; Norris, Marani, & Bonnell 
2000), with the soft band being from 25-50 keV and the hard band being from 100-300 keV. 
The energy bands for the various experiments are often fixed, yet choices can be made that 
are reasonably close to this standard (see Table 5). The individual lags calculated for each 
burst are given in Table 4. Some GRBs (those from Konus and SAX) do not have two 
channel light curves available, while others are too faint and hence have the uncertainty in 
the lag too large to be useful. 

The variability (V) of a burst is a quantitative measure of whether its light curve is 
spiky or smooth. A reasonable measure of this is to calculate the normalized variance of the 
observed light curve around a smoothed version of that light curve (Fenimore & Ramirez- 
Ruiz 2000). Unfortunately, this one-sentence description has a number of free parameters 
that can greatly change the calculated V. Fenimore & Ramirez-Ruiz (2000) and Reichart et 
al. (2001) experimented with a variety of choices. Schaefer, Deng, and Band (2001) selected 
effective choices as being those that produced the best correlation between V and Ti ag for 
112 BATSE bursts, finding that the formulation of Reichart et al. was poor and that the 
final choices of Fenimore & Ramirez- Ruiz were best. As part of this paper, I have tried a 
wide variety of options in the definition of V and have settled for the definition which yields 
the least scatter in the correlation between V and the burst luminosity. In this definition, 
I take the smoothed light curve [C smooth to be the original light curve (C) smoothed with 
box-smoothing where the width of the box is 30% of the time for which the light curve is 
brighter than 10% of its peak flux. The counts per time bin in the background-subtracted 
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light curve is C with an uncertainty of ac, while the counts in the smoothed light curve are 
C smooth- With this, the variability is 

V=([(C- C smooth ) 2 - (y 2 c\/C 2 smooth , max ), (9) 

where C smoot h,max is the peak value of C smoot h, and the angle-brackets indicate an average 
over all time bins where the smoothed light curve is greater than 10% of C smoo th,max- For 
each satellite experiment, I have chosen a fairly broad energy band over which to construct 
the light curve (see Table 5). With this definition, I have calculated V for most of the GRBs 
and placed these values in Table 4. Again, some GRBs are not included because the light 
curve might not be available, the light curve might have significant time gaps, or the burst 
brightness might be so faint that the derived uncertainty on V was too large to be useful. 

The minimum rise time {trt) in the GRB light curve is taken to be the shortest time 
over which the light curve rises by half the peak flux of the pulse. For bright bursts, it is 
an easy calculation to search time intervals before each peak for the shortest one in which 
the rise is half the peak brightness. In cases where the rise from one time bin to the next is 
greater than half-peak, I take the rise time to be the appropriate fraction of the bin width. 
But for faint bursts where the normal background noise provides large Poisson fluctuations, 
this procedure will always find fast rise times, so as a remedy I always bin the light curve 
until the uncertainty in individual points in the light curve are 10% of the peak flux (or less). 
In this binned light curve the fluctuations to mimic a rise by 50% of the peak flux would 
have to be roughly 5-sigma, and this eliminates the effects of Poisson noise. A cost of this 
binning is that rise times much shorter than the binned time intervals cannot be measured. 
The energy bands used for these light curves are the same as used for calculating V (see 
Table 5). These rise times will depend on the exact choice of the first bin, so I have adopted 
the average of the derived rise times over all possible start bins as being the minimum rise 
time. The uncertainty in the minimum rise time is then the RMS of the values for all the 
start bins. These values are tabulated in Table 4. 

The number of peaks (N peak ) in a GRB light curve varies from one for a simple fast-rise- 
exponential-decay (FRED) light curve shape to a dozen or more for complex spiky bursts. 
Let the overall maximum in the background subtracted light curve to be C max . I take a peak 
to be any local maximum which rises higher than C ma:E /4 and which is separated from all 
other peaks by a local minimum that is at least C ma x/^ below the lower peak. Again, this 
is simple to calculate for bright bursts but is not realistic for faint bursts with substantial 
Poisson noise. This is solved by binning the GRB light curve such that each bin has a one- 
sigma uncertainty that is less than C max /12. As before, the N pea k value can change with the 
phasing of the binning, but here the number of peak usually does not change. The energy 
bands and time resolutions for the light curves are specified in Table 5. The resultant N peak 
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values are listed in Table 4. 

The luminosity relations will be power laws of either L or E 7 as a function of Ti ag , V, 
Ep e ak, trt, and N peak . Both L and E 1 will have to be recalculated with luminosity distances 
appropriate for every cosmology. Table 4 has compiled all the values needed for producing 
a GRB HD for any particular cosmology. Columns 1 and 2 are the GRB identifier and the 
redshift. Column 3 is Pboio (the flux over the peak 1 second interval from 1-10,000 keV in 
the rest frame of the GRB) and its uncertainty in units of erg/cm 2 /s. Column 4 gives S^zo 
and its uncertainty in units of erg/cm 2 . Column 5 lists the beaming factor, F^ eam . Column 

6 tabulates the lag time for each burst in units of seconds in the Earth rest frame. Column 

7 gives the V value. Column 8 lists the observed E peak value from Table 2, provided that a 
definite measure of E peak is known. Columns 9 and 10 present the values for trt in seconds 
in the Earth rest frame and N peak . 



3. Outliers 

In any large set of observations taken from widely heterogeneous sources, there will be 
outliers. Many of these outliers arise because of some sort of error, rather than simply being 
a tail of some ordinary distribution of measurement errors. If these errors are allowed to 
remain, then they will dominate the fits and bias the results. So the problem is to recognize 
the outliers without throwing out the error-free values. I will adopt the standard solution 
of throwing out luminosity indicators or GRBs only if they deviate from best fit (with the 
outlier included) by 3-sigma or more. For some of the outliers, I have rejected the use of the 
entire burst, while for others I will impeach only a specific measurement. Here is a detailed 
list of the outliers that I have rejected: 

GRB 980425 is the smooth and faint burst which was identified with the supernova 
1998bw in a nearby galaxy at z=0.0085. From measures for the lag, E peak) and the rise 
time, I get a combined distance modulus of 40.35 mag whereas the real distance modulus 
(independent of any choices for the cosmology) is 32.8 mag or so. As with previous workers, 
I find that GRB 980425 is a distant outlier. The reason is likely that this very low luminosity 
event is somehow greatly different than classical GRBs and so the luminosity indicators do 
not apply. 

GRB 990123 was the original optical transient observed during the burst by ROTSE 
(Akerlof et al. 1999). Its redshift is a confident 1.61, and this is in good agreement with 
all five luminosity indicators except for the rise time. The minimum rise time in this well 
measured light curve (2.3 s) is much too long for a burst of this high luminosity. In the 
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calibration curve for the Ti ag — L relation, GRB 090123 has its rise time as a distant outlier 
and this rise time is rejected. I can only suggest that some bursts will happen to have their 
collisions all at far past the minimum radius and hence will have a long rise time for their 
luminosity. 

GRB 020819 was optically dark yet had a long fading radio transient that provides an 
accurate position (Jakobsson et al. 2005b). The indicated position is on top of a "blob" 
which appears near but definitely outside a barred spiral galaxy measured to be at z = 0.41. 
No redshift for the "blob" has been measured. GRB 020819 is an outlier on the Ti ag , V, tr T , 
and N pea k calibration curves. The combined value is 44.27 mag while the barred spiral 
galaxy has a distance modulus of 41.74 mag. Based on the luminosity indicators, the GRB 
is likely at z > 1 with a value around 1.7 being preferred. I suggest that the "blob" is the 
real host galaxy at z ~ 1.7, and this prediction can be tested with a spectrum of the "blob". 

GRB 030328 has the same situation as GRB 990123, in that all five luminosity indicators 
are in good agreement with the observed redshift except for the rise time (3.9 ±0.2 s) despite 
the burst being bright enough that faster rises would be obvious in the light curve. The two 
GRBs even have similar light curves. Again, I will reject the rise time as being a distant 
outlier. 

GRB 031203 does not have an available light curve or an E peak value, so it cannot be 
included in my analysis. Nevertheless, this burst is an outlier for both the T\ ag — L and 
Epeak — L relations (Sazanov, Lutovinov, & Sunyaev 2004). The trouble is that the reported 
measure of Ti ag = 0.24 ± 0.12s and E peak > 190 keV are both signs of a high luminosity event, 
whereas the burst is low luminosity at z = 0.105. It cannot be that the redshift is greatly 
in error since an underlying supernova (SN 2003 ) has been measured spectroscopically 
(Tagliaferri et al. 2004). The general thinking is that GRB 031203 is another unusual low- 
luminosity event like GRB 980425 (Sazanov, Lutovinov, & Sunyaev 2004; Soderberg et al. 
2004; Watson et al. 2004) for which the luminosity relations do not apply. 

GRB 050315 is a far outlier for the E pea k — L and E pea k — E 1 relations, while it is on 
the edge for the trt — L and on the best fit calibration line for the V — L relation. The 
combined jj, is 43.71 mag while the distance modulus of z — 1.948 is 45.87 mag. The problem 
is unlikely to be an error in redshift. This is because the outliers must be resolved by shifting 
to lower redshift yet the measured z comes from absorption lines in the afterglow (Kelson & 
Berger 2005) and thus it cannot be lowered. A real problem for the E peak — E y relation is 
that the x-ray afterglow light curve has at least four breaks in it (Vaughan et al. 2006) so 
any assignment of tj et can only be problematic. Another big problem for both the E pea k — L 
and E pea k — E y relations is that Vaughan et al. (2006) reports the E pea k value (I6I7 8 keV 
with error bars at the 68% confidence level) to be at the edge of the spectral energy range 
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for Swift (15-150 keV). In this case, the detection of a break can only be problematic. Both 
the E peak — L and E peak — E 7 relations would be easily satisfied if E peak is from 80-150 keV 
at the other edge of the Swift spectrum. So my suggested solution is simply that E peak is 
really at ~ 150 keV where Swift cannot see it. Despite the likelihood of this solution, GRB 
050315 must remain a rejected outlier for my analysis. 

In all, I have rejected only three bursts as outliers plus two rise times as outliers. This 
is fairly good given the many bursts and observations going into Tables 2-4. 

4. Luminosity Relations 

The luminosity relations are connections between measurable parameters of the light 
curves and/or spectra with the GRB luminosity. Specifically, I will be using the power law 
relationships between r lag - L, V - L, E peak - L, E peak - E 7 , t rt - L, and N peak - L. This 
section will discuss the calibration of all six relations. 

The calibration will essentially be a fit on a log-log plot of the luminosity indicator 
versus the luminosity. For this calibration process, the burst's luminosity distance must be 
known to convert Pboio to L (or Si> i to E y ) and this is known only for bursts with measured 
redshifts. However, an important point is that the conversion from the observed redshift to a 
luminosity distance requires some adopted cosmology. This means that every cosmology will 
have a separate calibration. Fortunately, as we will see below, the calibration only weakly 
depends on the input cosmology. If we are interested in calibration for purposes of GRB 
physics, then it will be fine to adopt the calibration from some fiducial cosmology (say, the 
concordance cosmology with w = —1). But if we are interested in testing the cosmology, then 
we have to use the calibration for each individual cosmology being tested. For a particular 
cosmology, the theoretical shape of a HD has to be compared with the observed shape when 
the burst distances are calculated based on calibrations for that particular cosmology. In 
practice, this means that the model and observed luminosity distances are compared in a 
chi-square sense as cosmological parameters vary with the observed values changing with the 
cosmology. Thus, any test of cosmological models with a GRB HD will be a simultaneous 
fit of the parameters in the calibration curves and the cosmology. 

A comparison of this case with that of the supernova HD is instructive. Schematically, 
the calibration of the luminosity relation for supernovae (i.e., decline rate versus luminosity) 
can be accomplished for nearby events for which the luminosity distances have no dependance 
on the cosmological parameters. Once this calibration is learned from nearby events, then 
it can be applied to distant supernovae with confidence. In this case, the calibration does 
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not depend on cosmology and the deduced positions of the high redshift supernovae in a 
Hubble diagram will not depend on the cosmology being tested. But in practice, many of 
the supernovae in Hamuy et al. (1996) are sufficiently distant that small effects of varying 
cosmology must be introduced, and then even the supernova HD becomes a simultaneous 
fit involving both calibrations and cosmologies. Another feature of the supernova HD case 
is that it involves both nearby events where the distances come from Cepheids that are 
independent of the Hubble constant (H ) and far events where the distances depend on the 
Hubble constant. In such a case, distortions can appear in the shape of the HD depending on 
the adopted value for the Hubble constant. To solve this problem, Perlmutter et al. (1997; 
1999) adopted a method where the " Hubble-constant-free luminosity distance" is used along 
with the calibration with a " Hubble-constant-free B-band absolute magnitude at maximum" , 
where they subtract 5 log H from the peak absolute magnitude (M) and add it back in again 
into the distance modulus (//) in the usual equation for the observed magnitude (m). Thus, 
m = M + fi becomes m = (M — 5 logiJo) + (// + 5 logifo)- But with GRBs, we do not have 
the problem of mixing distances that are both dependent and independent of the Hubble 
constant. So for GRBs, the "Hubble-constant-free" equation gives identical results as the 
normal equation. 

The conversion from redshift (z) to luminosity distance (di) depends on the cosmology. 
Throughout this paper, I will assume that the Universe is flat, as indicated by strong sets of 
theoretical and observational arguments. For a flat Universe with the concordance cosmology 
of w — — 1, this relation is 



Here, VLm is the dimensionless matter density and Vt\ = 1 — VLm for a flat Universe. If 
we allow for the possibility that the Dark Energy changes with time (i.e., w varies with 
redshift), then we get a different relation. In the absence of any real understanding of how 
Dark Energy changes (c.f. Szydlowski, Kurek, & Krawiec 2006), a reasonable approach 
for analysis of observations is to simply expand w as some function of redshift. Riess et al. 
(2004) introduced the expansion w(z) = wo+w'z, but this has been widely criticized as being 
unsuitable at high redshifts because an exponential term grows large for z > 1. A better 
expansion (i.e., one that is well bounded at high redshift) is to adopt w(z) = WQ+Wazil+z)^ 1 
(Linder 2003). This gives 



d L = cH \l + z) / dz'[(l + z') 3 Q M + Oy{(1 + ^) 3 ( 1 +-o+- a )} e -3« )a2 /(i+,) ] -i/2 > (n) 



Jo 

All such expansions are attempts to parameterize the redshift dependance of w with a min- 
imum of adjustable parameters. 




(10) 
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In addressing the question of testing whether the Dark Energy changes with redshift, it 
will be convenient to compare two cases: one representing an unchanging w — — 1, while the 
other is representative of some reasonable case of a variable w. Such a comparison can be 
helpful for seeing how much things change (e.g., the relative positions of the observed GRBs 
in the HD) between the constant and variable cases. The only 'representative' case for a 
variable w is that of the best fit for the 'gold sample' of supernovae (Riess et al. 2004) with 
w = —1.31 and w' = 1.48. In this paper, I will only use the w(z) = w + w'z expansion for 
purposes of making this comparison. For all other purposes (like constraining the possible 
variations in the Dark Energy in Section 7), I will adopt the w(z) = wo + w a z{l + z)^ 1 
expansion of Linder (2003). 

The observed luminosity indicators will have different values from those that would 
be observed in the rest frame of the GRB. That is, the light curves and spectra seen by 
Earth-orbiting satellites suffer time-dilation and redshift. The physical connection between 
the indicators and the luminosity is in the GRB rest frame, so we must take our observed 
indicators and correct them to the rest frame of the GRB. For the two times (ri ag and trt), 
the observed quantities must be divided by 1 + z to correct for time dilation. The observed V 
value varies as the inverse of the time stretching, so our measured value must be multiplied 
by 1 + z to correct to the GRB rest frame. The observed E peak value must be multiplied by 
1 + z to correct for the redshift of the spectrum. The number of peaks in the light curve is 
defined in such a way as to have no z dependance. The dilation and redshift effects on 9j et 
and E 1;iso have already been corrected in equations 1 and 2. A possibly substantial problem 
for the Ti a g, V, and t rt relations is that we are in practice limited to the available energy 
bands (c.f. Table 5) whereas these correspond to different energy bands in the GRB reference 
frame. Ideally, we would want to measure these indicators in observed energy bands that 
correspond to some consistent band in the GRB frame. 

The luminosity relations are all expressed as power laws and thus should be a linear 
function of the logarithms of the corrected indicator versus the logarithms of the burst 
luminosity. These are displayed in Figures 1-6 for the concordance cosmology (w = — 1 and 
Q = 0.27 in a flat universe). The scatter of the data about the best fits is consistent with a 
Gaussian distribution in this log-space. 

The plots of the luminosity relations show significant error bars in both the horizontal 
and vertical axes. We would get different best fit slopes depending on whether we calculate 
a standard ordinary linear regression where the residuals in the luminosities are minimized 
or a linear regression where the residuals in the luminosity indicators are minimized. The 
measurement uncertainties in both the luminosities and their indicators is smaller than the 
observed scatter about the best fit lines. This means that there must be some additional 
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source of intrinsic scatter which is dominating over the simple measurement errors. And 
this scatter appears to be independent of both luminosity and redshift. In this case, the 
conclusion is to perform the ordinary least squares without any weighting. (The use of 
weighted least squares results in almost identical best fits.) The best fit values (from ordinary 
least squares with no weighting) are given in many reference books (e.g., Press et al. 1992). 
In the luminosity relations, the two variables are not directly causative (with both being 
determined by Tj et ), so the bisector of the two ordinary least squares (Isobe et al. 1990) 
will be used. To be specific with the lag-luminosity relation as an example, with Yi = log L 
and Xi = \og(r lag ) — log(l + z), we are seeking the best fit to Yi — a + bXi. The simple 
averages of Xi and Y { are denoted as X and Y. The intercept is a = Y — bX. For 1^0, 
the best estimate of a will be sensitive to small errors in b and the confidence region for a 
and b will be a long thin ellipse. To avoid these problems, we can recast the equation to be 
fit as Yi — Y — b(Xi — X). 

In the following sections, I will present each of the luminosity relations and their derived 
best fit equations. Also, for each relation, I will sketch out the theoretical justification 
and the physics behind each relation. (In fact, three of these relations were theoretically 
predicted before they were confirmed with data.) Despite have a theoretical understanding, 
the relations are essentially empirical, with the tightness of the relation being more important 
than any physical model. This case is also true for the calibration of Type la supernovae. 

4.1. Lag versus Luminosity 

The GRB lag time was first identified as a luminosity indicator by Norris, Marani, 
& Bonnell (2000) based on six BATSE GRBs with optical redshifts. This relation has a 
closely inverse relation between luminosity and lag, with very long lag events being very 
low luminosity and very short lag events being very high luminosity. This relation (and the 
V — L relation) was confirmed by Schaefer, Deng, & Band (2001) based on the predicted 
T iag ~ V relation appearing for 112 BATSE bursts with no measured redshifts. Norris (2002) 
demonstrated that the long lag bursts (hence the lowest luminosity and generally the closest) 
show a concentration to the supergalactic plane. 

The physics of the lag-luminosity relation is apparently simple (Schaefer 2004), and 
indeed the relation should have been deduced a decade ago as a necessary consequence of 
the empirical/theoretical Liang-Kargatis relation (Liang & Kargatis 1996; Crider et al. 1999; 
Ryde & Svensson 2000; 2002). The Liang-Kargatis relation is that the time derivative of 
Epeak is proportional to luminosity, as established by empirical measures or for the expected 
case where the cooling of the shocked material is dominated by radiative cooling. Schaefer 
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(2004) demonstrated that this time derivative equals a constant divided by the lag, so that 
we deduce L oc r z ~*. In essence, the shocked material will cool off at a rate dictated by its 
luminosity. The lag is related to the time for a pulse to cool somewhat; if the burst is highly 
luminous then the cooling time (and hence the lag) will be short, while if the burst has a 
low luminosity then the cooling time and the lag will be long. 

The calibration data is plotted in Figure 1 as log[r/ aff /(l + z)\ versus L. This plot is 
based on an assumed concordance cosmology (w = —1), JIm — 0.27 (and hence Q\ = 0.73), 
and equation 10. The best fit linear regression line is also plotted. The equation for this 
calibration line is 

logL = 52.26 - 1.01 logfa ag (l + z)-70.1s]. (12) 

I note that the slope in this relation is satisfactorily close to the theoretical value of -1. The 
one-sigma uncertainties in the intercept (a = 52.26) and slope (b = —1.01) are a a = 0.06 
and ay, = 0.05. 

What is the uncertainty on any L value deduced from this relation? There will be two 
components that must be added in quadrature. The first is the normal measurement error 
obtained by propagating the uncertainty through the last equation. The second will be the 
addition of some systematic error that accounts for the extra scatter observed in Figure 1. I 
will label this systematic error as <Ji a g,sys- Then, we have 

<l = °l + (?b logfr as (l + *)-70.1s]) 2 + (0AM3ba lag /T lag ) 2 + <j 2 lag , sys . (13) 

The value of (Jiag,sys can be estimated by finding the value such that a chi-square fit to the 
lag-luminosity calibration curve produces a reduced chi-square of unity. With this, I find 

&lag,sys 0.39. 

How much does this calibration depend on the cosmology? To get an idea of the char- 
acteristic changes, we can compare the calibration for the concordance w — — 1 cosmology 
(i.e., equation 10) with that derived for the changing dark energy case given as the best fit to 
the 'gold sample' of supernovae (w = —1.31 and w' = 1.48) from Riess et al. (2004). This 
best fit cosmology from Riess produces calibration intercept and slope of 52.18 and -0.96. As 
such, we see that the variation in the calibration of the lag-luminosity relation only weakly 
depends on typical range of cosmological models. 

4.2. Variability versus Luminosity 

The GRB variability (V) was first recognized as a luminosity indicator by Fenimore & 
Ramirez-Ruiz (2000) based on seven BATSE GRBs with optical redshifts. The validity of 
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the V — L relation was confirmed with additional bursts by Reichart et al. (2001) and by the 
existence of the predicted T\ ag — V relation seen with 112 BATSE bursts without redshifts 
(Schaefer, Deng, & Band 2001) and as updated to 551 BATSE bursts without redshifts 
(Guidorzi 2005). Lloyd- Ronning & Ramirez- Ruiz (2002) also confirmed the existence of the 
V — L relation by demonstrating the existence of the predicted V — E peak relation for 159 
BATSE GRBs with no known redshift plus 8 BATSE GRBs with optical redshifts. Fenimore 
& Ramirez-Ruiz (2000), Reichart et al. (2001), Schaefer, Deng, k Band (2001), Guidorzi et 
al. (2005), and Li & Paczyhski (2006) have presented a variety of definitions of V which use 
various smoothing functions and parameters and choice of time intervals. 

The origin of the V — L relation is based in the physics of the relativistically shocked 
jets. Detailed explanations within this model have been provided by Meszaros et al. (2002) 
and Kobayashi, Ryde, & MacFadyen (2002). In general, both V and L are functions of the 
bulk Lorenz factor of the jet (Tj et ), where the luminosity varies as a moderately high power 
of Yj et and where the high values of Tj et make for smaller regions of visible emission and 
hence the faster rise times and shorter pulse durations result in high variability. 

The calibration plot for the V — L relation is given in Figure 2 along with the best fit 
line. This best fit can be represented with the equation 

logL = 52.49 + 1.771og[F(l + z)/0.02]. (14) 

The one-sigma uncertainties in the intercept and slope are a a = 0.22 and a b = 0.12. The 
uncertainty in the log of the luminosity is 

<l = °l + (<r b log[V{l + z)/0.02]) 2 + (0A343ba v /V) 2 + a^ sys . (15) 

The chi-square of the points about the best fit line is unity when av, sys = 0.40. 

The best fit cosmology from Riess et al. (2004) produces calibration intercept and slope 
of 52.22 and 1.67 respectively. This alternate cosmology changes the calibration parameters 
by about one-sigma. 



4.3. E peak versus Luminosity 

E peak has been strongly correlated with both L (Schaefer 2003b) and E ltiso (Amati et 
al. 2002). The two relations are likely caused by different physics, with the E peak — E ltiso 
relation being an approximation of the E peak — E 1 relation (discussed in the next subsection) 
that is related to the total energetics of the burst. The E peak — L relation is different because 
it is related to the instantaneous physics at the time of the peak. The idea is that the peak 
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luminosity varies as some power of Tj et while the E peak also varies as some other power of Tj et , 
so that E peak and L will be correlated to each other through their dependance on Tj et . Indeed, 
a detailed analysis of the situation resulted in a prediction that L oc [E peak (l + z)] N ^ M+1 ^ 
with N/(M + 1) ~ 2.5, and this prediction was confirmed (Schaefer 2003b) from sets of 20 
and 84 bursts. This analysis also explains why the observed distribution of E peak can possibly 
be so narrow (Mallozzi et al. 1995) despite very wide ranges in both Tj et and 1 + z, as well 
as explains why the observed average E peak varies as a particular function of the observed 
peak flux (Mallozzi et al. 1995). 

Figure 3 plots the values of log[E peak (l + z)\ versus logL for all bursts with available 
data. The best fit line is plotted and is represented by the equation 

logL = 52.21 + 1.681og[-E7peafc(l + z)/ZQIQkeV\. (16) 

The one-sigma uncertainties in the intercept and slope are a a = 0.13 and = 0.05. The 
uncertainty in the log of the luminosity is 

<l = °l + {<J b \og[E peak {l + z)/?>mkeV}) 2 + (0A343ba Epe jE peak f + a\^ sys . (17) 

The chi-square of the points about the best fit line is unity when o~E peaktS ys — 0.36. 

With Riess' best fit cosmology (wo = —1.31 and w' = 1.48), the intercept and slope are 
52.11 and 1.60 respectively. Again, the calibration is only weakly dependent on the input 
cosmology. 

4.4. E peak versus E 1 

Ghirlanda, Ghisellini, & Lazzati (2004) discovered a tight correlation between E peak and 
E 1 . This is an improvement on (and combination of) both the E 1 = constant relation of 
Bloom et al. (2003) and the E peak — -E 7i j SO relation of Amati et al. (2002). The physics of 
the E peak and E 1 relation is well explained as a simple consequence of viewing geometry and 
relativistic effects within a standard jet model (Eichler & Levinson 2004; Yamazaki, Ioka, 
& Nakamura 2004; Rees & Meszaros 2005; Levinson & Eichler 2005). This relation has the 
advantage of being one of the tightest for GRBs. But this relation can only be used for the 
minority of GRBs with redshifts because a jet break has to be identified and measured in the 
afterglow light curve amongst the many bumps and breaks that are ubiquitous in afterglows. 

Figure 4 shows the calibration curve for the E peak — E^ relation. The best fit is 



log£ 7 = 50.57+ 1.63 \og\E peak (l + z)/300ifeeV]. 



(18) 
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The one-sigma uncertainties in the intercept and slope are a a = 0.09 and a b = 0.03. The 
uncertainty in the log of the burst energy is 

<e 7 = °l + (a b \og[E peak (l + z)/300keV]) 2 + (0A3A3ba Epe jE peak ) 2 + a\^ ys . (19) 

The chi -square of the points about the best fit line is unity when &E 1:S ys — 0.16. 

With the Riess best fit cosmology, the intercept and slope are 50.50 and 1.59 respectively. 

4.5. t rt versus Luminosity 

The variability of a light curve is a peculiar construction as we have no clear idea of 
what we are trying to measure and it is difficult to understand the physics of 'variability'. In 
an effort to understand the meaning of variability, I calculated variability for a wide range of 
simulated light curves constructed from individual pulses with shapes as given by the model 
in Norris et al. (1996). The most important determinant of the V value was the rise time in 
the light curves, with other properties (like fall time, pulse duration, burst duration) being 
of lesser importance. Thus, it seems that the V might be just a surrogate measure for the 
rise time in the light curve. And the rise time can be directly connected to the physics of 
the shocked jet. Indeed, for a sudden collision of a material within a jet (with the shock 
creating an individual pulse in the GRB light curve), the rise time will be determined as the 
maximum delay between the arrival time of photons from the center of the visible region 
versus the arrival time of photons from the edge of the visible region. Because the jet material 
is traveling at very close to light speed, this delay time is simply from the longer path length 
traveled (much like an echo). This delay depends on the angular size (as viewed from the 
center of the GRB) of the visible region, which then depends on Tj et , so that rise times 
during a burst are proportional to Tj^ t . The proportionality constant depends on the radius 
from the GRB that the material is shocked. But there is some minimum radius under which 
material is optically thick and inefficient at radiating and hence faint. This minimum radius 
is roughly a constant from burst-to-burst (Panaitescu & Kumar 2002), and this means that 
the minimum rise time should be roughly proportional to 17^. The scatter in this relation 
will depend on how close the collisions in the jet occur to the minimum radius, so we could 
expect a substantial amount of scatter. The burst luminosity also scales as T^ et for N ~ 3. 
With both trt and L being functions of Tj et , I predicted that the minimum rise time should 
be a luminosity indicator with L oc r R ^ 2 (Schaefer 2002). 

This prediction has proven true, for example as shown in Figure 5. The best fit is 



logL = 52.54- 1.21 log[T fiT (l + z^/OAs}. 



(20) 
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The one-sigma uncertainties in the intercept and slope are a a = 0.06 and a b = 0.06. The 
uncertainty in the log of the burst energy is 

<l = °l + K log[rKr(l + zy'/O.ls}) 2 + (0A343ba TRT /r RT ) 2 + a 2 RTjSys . (21) 

The chi-square of the points about the best fit line is unity when <jRT, sy s = 0.47. 

With the Riess cosmology, the intercept and slope are 52.42 and -1.14 respectively. 

4.6. N peak versus Luminosity 

The number of peaks in a light curve depends on how many collisions between packets 
of material in the jet occur during the duration of the burst. This number will be deter- 
mined by many factors, including the exact realization of turbulence in the source and the 
distribution of velocities and densities in the jet. However, some of the individual peaks 
might occur sufficiently close in time that these peaks will appear as one. If the individual 
pulse durations are somewhat longer than the separation in time, then the two pulses will 
not be distinguishable as being separate. The pulse durations (D pu i se ) scale as the rise times 
(Nemiroff 2000) and hence will scale as Tj 2 t or L~ 2 / N (see previous subsection). Schaefer 
(2003b) has presented several theoretical and observational arguments with the combined 
result that N = 3.14 ± 0.34. For high luminosity bursts all collisions will result in distinct 
pulses in the light curve, while low luminosity events will have many of the collisions resulting 
in overlapping broad pulses. Thus, a burst with many peaks can only be a high luminosity 
event because this is the only way to get narrow peaks that avoid merging together. A 
burst with one or a few peaks could either be high luminosity (with few shell collisions) or 
low luminosity (with all the collisions producing merged peaks). The maximum number of 
distinct peaks in a light curve with a duration of T go is 1 + (T go / D pu i se )e~^ 1+Dpulse ^ Tgo ^ for 
Poisson distribution of collisions. With D pu i se /T 90 ~ [Lj L )~ 2 / N and N ~ 3.14, we can 
translate the observed number of peaks in a burst (N pea h) into a lower limit on the burst 
luminosity. Again, this analysis was a theoretical prediction that was tested and shown to 
be true (Schaefer 2002). 

Figure 6 has a plot of the burst luminosity versus the number of peaks in a burst's light 
curve. The curved line is the theoretical limit, and this was drawn for logLo = 50.08 so as 
to tuck the limit up against the lower envelope of the data. The important point is that 
there is an observational lower limit that corresponds in shape to the theoretical prediction, 
and this confirms the general picture. 

The theoretical limit is actually quite straight when plotted in log-log space. So a 
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convenient representation would be a simple power law 

log L > 50.32 + 2 log[N peak ] for N peak > 2. (22) 

For Npeak — 1, there is no lower limit on the luminosity. This approximation is plotted in 
Figure 6 as the thin straight line segment. 

As a limit, equation 22 has little utility for Hubble diagram purposes. This means that 
I will only be using the first five luminosity indicators (sections 4.1-4.5) for the rest of this 
paper. However, there is utility in the N pea k — L relation for two other purposes. First, 
it provides a startling theoretical prediction and observational correlation that should be 
studied for further insights into burst physics. Second, counting the number of peaks is a 
fast way to spot some high luminosity GRBs. For example, if Swift sees a faint burst with 
many peaks, then the burst must be at high redshift. 



4.7. Combining the Derived Distance Moduli 

For every L or E 7 that we calculate, we can also derive a luminosity distance from the 
inverse-square law. The equations for this are 

d L = [L/(4nP Mo )] 1/2 , (23) 
d L = [£ 7 (1 + z)/{AixF beam S bolo )} 1 / 2 . (24) 
A distance modulus (//) can be calculated for every estimated luminosity distance as 

^ = 51og(d L )-5, (25) 

with di expressed in units of parsecs. The propagated uncertainties will depend on whether 
Pboio or Sboio is used; 

a, = [(2.5<7 logL ) 2 + (1.086a Pbo J P^) 2 } 1 ' 2 , (26) 

or 

<t„ = [(2.5a log£7 ) 2 + (1.086^/S^) 2 + {imQa Fbaa jF beam ) 2 ] 1/2 . (27) 
With five luminosity indicators, each burst will have up to five measured distance moduli 
and their one-sigma uncertainties, which I will label as fi± ± <r m , ^ ± °7i 2 > ^3 ^ °>3 > ^4 ^ °W > 
and /i 5 ± <r^ 5 for the five indicators from subsections 4.1-4.5 respectively. 

The best estimate \x for each GRB will be the weighted average of all available distance 
moduli. Thus, the derived distance modulus for each burst will be 

/^ = E^/<]/E<]' ( 28 ) 

i i 
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and its uncertainty will be 

^ = E<r 1/2 > (29) 

i 

where the summations run from 1-5 over the indicators with available data. The weighted 
average formalism takes care of the case where the various indicators have greatly different 
scatter. For example, if a GRB only has a distance modulus from the V — L and E pea k — E 1 
relations (i.e., yU 2 and /i 4 ), then the disparate uncertainties (with <t M2 <7 M4 ) will result in 
the variability contributing little to the final answer. 

A potential problem with equation 29 arises if the luminosity relations are correlated. 
In an extreme case, where two relations are perfectly correlated, the use of equation 29 
would incorrectly double the weight of what is really just one relation, which is equivalent to 
reducing the error bars by a factor of 1.4. This would be non-optimal but not disastrous. A 
less extreme example of this would be if the E peak — E lAso relation of Amati et al. (2002) was 
used at the same time as the E pea k — E 1 relation, as much the same physics and input are 
used in both. For application to this paper, I wonder whether the E pea k — L and E pea k — E 1 
relations are significantly correlated as they share a common input of E pea ^. To test this, I 
have tabulated /ij — fi z for all bursts and relations, where /j, z is the distance modulus derived 
from the observed redshift and some fiducial cosmology (here, the concordance cosmology). I 
find that the E pea k — L and E pea k — E 1 relations have a near-zero correlation coefficient, which 
demonstrates that the two relations are not correlated and hence independent. However, 
there is one significant correlation, with the correlation coefficient equalling 0.53 between the 
V — L and r RT — L relations. This is not surprising, as the V of a light curve is dependent 
on the rise time, and indeed the minimum rise time was originally proposed as a measure of 
'variability' that has a direct physical interpretation. From the extent of the fi2 — fi z versus 
A*5 — scatter plot, the correlation is roughly half of the intrinsic scatter. For the 71% 
of the bursts with both V and trt, the statistical weight of the V — L relation should be 
roughly cut in half, which is equivalent to the combined a M2 , being 20% too small. These 
two relations provide relatively little weight for well-observed bursts (with all five luminosity 
indicators known), so that the a M will only be 4% too small in these cases. 

The statistical weight of each measured distance modulus is a^, so the total statistical 
weight for each relation contributing to the HD will be a summation of <r Mi over all bursts. 
This will give us an idea of the relative contribution of each of the relations. I find the 
summed weights to be 30.3, 27.0, 54.2, 64.9, and 38.4 for the Ti ag — L, V — L, E peak — L, 
E peak — E 1 , and t rt — L relations respectively. This corresponds to percentages of 14%, 13%, 
25%, 30%, and 18% respectively. The E peak — E 1 relation is the most accurate of the five, 
but it does not dominate since only 27 bursts are useable. The V — L, E pea j, — L, and trt — L 
relations all have information for almost all bursts, so their percentages is an indication of 
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the size of the derived uncertainties. In all, we see that all five relations have comparable 
total weight, with none dominating and none being negligible. 

Each of these measures carries information, so it would be wrong to not use them with 
appropriate error bars. Optimally, we want to use all information, so that the use of any one 
indicator is neglecting the majority of the constraints. For example, will improve from 
typically 0.7 mag if only the E peak — E 1 relation is used to around 0.4 mag if all five relations 
are used. 

To illustrate this calculation and some intermediary values, I present the various values 
for pL in Table 6. The first column is the GRB identifier and the second gives the redshift. 
Columns 3-7 are the //j ± values all for the concordance cosmology case with Qm = 0.27 
and w = —1. In column 8, I give the derived distance modulus from equations 25 and 26. For 
comparison, column 9 gives the derived distance modulus for the case of Riess' cosmology 
(wq = —1.31 and w' = 1.48). In comparing the last two columns, there are shifts in the 
derived ji values depending on the assumed cosmology, in this case with a typical shift of 
0.25 mag with only a small and loose dependence on redshift. 

5. Hubble Diagram 

The GRB Hubble diagram is a plot of the distance moduli versus the redshifts. Table 
6 has already listed ji and z values for two cosmologies. As such, Figures 7 and 8 are GRB 
HDs as taken straight from Table 6. 

Figure 7 shows the GRB HD for 69 GRBs out to redshift >6 for the concordance 
cosmology (indicated by the curved line). Several aspects of this Figure are striking: First, 
the GRBs define a well-behaved curve. Second, the regime covered by the supernovae (z < 1 
except for 10 events with 1 < z < 1.7) is only a fraction of the left hand side of this plot. 
Third, about half the points are at z > 2 with good coverage spanning a redshift range 
far higher than that available even with satellite observations of supernovae. Fourth, the 
observed data points are in good agreement with the model curve (with a reduced chi-square 
of 1.05). Fifth, the implication is that the GRB HD is consistent with w — — 1 and an 
unchanging Dark Energy. 

Figure 8 shows the HD for the same 69 GRBs except that the luminosity relations 
were calibrated with a particular cosmology that has the Dark Energy changing with time. 
The goal here is to illustrate the relative lack of shifting of the fi values as the cosmology 
parameters change over a typical range. Between Figures 7 and 8, the points shift by an 
average of 0.23 mag (but this constant shift does not affect the shape of the experimental HD 
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curve) while the RMS scatter in the shift is only 0.10 mag largely independent of redshift. 
These small relative shifts are to be compared to the change in the model \x by 0.65 mag at 
z — 4. The conclusion is that the position of GRBs in the HD is independent of the input 
cosmology (over a reasonable range of parameters) to first order. The particular cosmology 
in Figure 8 is that of the best fit by Riess et al. (2004) as based on their 'gold sample' 
of supernovae. Riess' cosmology is seen to lie significantly below the data for z > 3 and 
above the data for middle redshifts of 1 < z < 3. This situation can be quantified with a 
chi-square parameter comparing the observations and the model. Riess' cosmology gives a 
chi-square of 80.3 while the concordance cosmology gives a chi-square of 72.3. This difference 
of chi-square implies a nearly three sigma rejection of Riess' cosmology when compared to 
the concordance cosmology. 

An important point is that the luminosity relations have to be calibrated for every set 
of cosmological parameters under consideration. Schematically, this is different from the 
supernova HD where the calibration of the decline rate versus luminosity relation can be 
done at low redshift independent of cosmology and then applied to high redshift events. 
All GRB distances are dependent on the cosmology because there is no low-redshift set to 
perform the cosmology-independent calibration. This forces us to make a simultaneous fit 
of the cosmology and the luminosity relations (Schaefer 2003a). This simultaneous fitting 
allows us to completely avoid any circularity. 

The procedure of calibrating the luminosity relation for each cosmology ensures that 
the constructed HD has the deviations between observation and model being near zero. 
The vertical positioning of the GRB points in the HD are thus 'normalized' out simply by 
increasing or decreasing the adopted luminosities. But this is fine, since it is the shape of 
the HD that provides the cosmology information. As the cosmological parameters vary, the 
deduced d L for each burst will vary, the positions in the calibration plots (Figures 1-5) will 
shift, the fitted calibration equations (equations 12, 14, 16, 18, and 20) will change, and 
the derived distance modulus in the HD will shift. In the change from the concordance 
to the 'Riess' cosmology, the model curve has moved downward with a shift that increases 
to higher redshift. As such, the value for each burst is lower for the 'Riess' cosmology 
and so the intercept in all the luminosity relations will decrease. An overall lowering of the 
intercepts will not change the shape of the HD. So what matters are the shifts relative to this 
overall offset. Fortunately, there is a strong mixture of burst luminosities and redshifts, as 
this avoids a degeneracy where the luminosity relations and the cosmology cannot become 
untangled. The fact that the highest redshift GRBs can only be high luminosity implies that 
there will be some correlations. The typical degree of this correlation has already been seen 
in the comparison of the \x values between the concordance and 'Riess' cosmologies. Here, the 
average scatter in the shifts (0.10 mag) is a small fraction of the shift in the model (0.65 mag 
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at z — 4), while the systematic shifts with z are also relatively small (0.1 mag from redshifts 
1 to 4). Another way to show that there is a good mixing between redshift and luminosity is 
that the median redshifts in the ranges 51 < logL < 52, 52 < logL < 53, and 53 < logL are 
1.57, 2.42, and 1.48 respectively. It is this mixing of redshifts and luminosities that allows 
for the luminosity relations to distinguished from the cosmology during a simultaneous fit. 

What is the accuracy of a single GRB in the HD? One measure of this is the median 
a^, which is 0.60 mag. (The sizes of values are apparently reliable since the reduced 
chi-square for the HD is near unity.) Another measure is that the deviations have an RMS 
scatter of 0.65 mag. 

A measure of the accuracy in the GRB HD is to plot my derived redshifts (based on the 
five luminosity indicators) versus the known redshifts (based on optical spectroscopy). This 
is shown in Figure 9. The scatter of the derived redshifts about the diagonal line (where 
the derived redshifts equal the known redshifts) has a reduced chi-square near unity, and 
this means that the quoted error bars are realistic. If we look at the differences between the 
derived and known redshifts, the RMS error is 26%. The maximum error is 69%, with the 
second worst error being 56%. One the good side, 15% of the Swift bursts have 14%-18% 
fractional errors in redshift. 

The number of degrees of freedom in a chi-square comparison of the GRB data and 
theoretical models depends on the question being asked. If we want to do an evaluation of 
the concordance cosmology, then there are no free parameters for fitting. (The parameters 
for fitting the luminosity functions do not count as we are not optimizing on the goodness of 
fit in the HD. If we instead let the calibration parameters vary so as to optimize the fit in the 
HD, then we get correspondingly lower chi-squares in the HD yet poor calibration fits.) That 
is, we just accept Qm — 0.27, Q\ = 0.73, and w = —1 and we have no changeable parameter 
to optimize the goodness of fit in the HD. In this case, the number of degrees of freedom is 
just the number of bursts in the GRB HD (69). If we are seeking to constrain cosmological 
parameters, then the number of degrees of freedom will be the number of bursts (69) minus 
the number of adjustable parameters. So if we take a flat universe with fl M = 0.27 as a prior 
while asking for the best values of w and w a , then we have 69-2=67 degrees of freedom. 
If we test a more general model where Qm, ^a, an d w a are allowed to vary so as to 
optimize the fit in the HD, then we have 65 degrees of freedom in the chi-square. 

Two further tasks need to be accomplished before the GRB HDs, such as presented 
in Figures 7 and 8, can be used to constrain cosmological parameters. First, we have to 
evaluate biases arising from gravitational lensing and Malmquist effects. Second, we have to 
run chi-square fits over ranges of parameters and marginalize 'nuisance' parameters. 
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6. Gravitational Lensing and Malmquist Biases 

The light from high redshift GRBs travels a long way in our Universe, and there is 
significant probability that it will pass close enough to a galaxy to be gravitationally lensed. 
Infrequently, this lensing will magnify the apparent brightness of the burst as seen here at 
Earth, perhaps with high amplification. More commonly, lensing will demagnify the image 
and make it appear somewhat fainter. The average magnification must be exactly unity. This 
magnification and demagnification can cause scatter in the observed HD. To see this, imagine 
some particular GRB at some particular redshift with some particular luminosity indicators. 
In the case with no lensing, the luminosity indicators would return the approximately correct 
luminosity and we would use the observed peak flux to return nearly the correct luminosity 
distance and \x for plotting on the HD. In the case with lensing, neither the observed redshift 
nor the luminosity indicators will be any different, but the observed peak flux will have 
changed by some factor. With a shifted peak flux, the deduced luminosity distance and \x 
will be shifted and the GRB will be plotted on the HD with a somewhat shifted position. 
Most bursts will appear slightly fainter than then they should if unlensed, but some will 
appear brighter and rarely they can appear much brighter. There is no observational means 
of measuring the magnification for any GRB. The result will be an increased and irreducible 
scatter in the HD, and indeed this expected scatter is likely to account for some of the 
systematic scatter observed in the calibration plots (Figures 1-5). 

Several groups have calculated the probability distributions for magnifications as a func- 
tion of redshift. Holz & Linder (2005) calculate that the effective RMS of the scatter from 
lensing is roughly 0.093,2 mag up to z = 2. They find that with 70 sources at z — 1.5 the 
lensing effects will average out to errors of less than one percent. Premadi & Martel (2005) 
present probability distributions for magnifications up to a redshift of 5, and find that the 
dispersion turns over for z > 1 so that the Holz & Linder result should not be used for z > 1. 
Oguri & Takahashi (2006) make a first attempt to calculate the effects of lensing on a GRB 
HD, and they conclude that the lensing effects are negligible. All three results point towards 
lensing not being a significant bias for the GRB HD. 

The effects on the HD are more complicated than just the addition of some scatter in the 
deduced fx. The reason is that there is effectively a threshold in burst apparent brightness for 
its being detected, its redshift being measured, and luminosity indicators being measured. 
Although it is a 'fuzzy' limit, there will be a corresponding distance limit beyond which GRBs 
will not be included in the sample in Table 2. With gravitational lensing, bursts just outside 
this threshold might be magnified in brightness and included in this sample, whereas bursts 
just inside this threshold might be demagnified in brightness and excluded from this sample. 
This threshold effect will result in the highest redshift events appearing systematically more 
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luminous than they would be in the absence of lensing. The quantitative evaluation of these 
effects depends on the shape of the luminosity function (Oguri & Takahashi 2006). 

There is a similar set of effects (referred to generically as Malmquist bias) that result 
from observational uncertainties and geometry effects near the threshold of detection. With 
observational errors and intrinsic scatter in the luminosity relations, bursts just inside the 
distance limit will be excluded if the random fluctuations push the apparent brightness below 
the detection threshold, while bursts just outside the distance limit will be included if the 
random fluctuations push the apparent brightness above the detection threshold. There is 
more volume just outside the distance limit than there is just inside, so this will result in 
more over-bright GRBs being included than under-bright GRBs being excluded. A similar 
effect occurs for the luminosity because there are always many more lower-luminosity events 
that can have random fluctuations excluded in the sample than there are higher-luminosity 
events that can have random fluctuations included in the sample. 

Gonzalez & Faber (1997) provide a general expression that can be used to calculate 
both the gravitational lensing biases and the Malmquist biases. In general for any particular 
GRB, the expression for the bias is 

fl + B)D ~ (r\D) - $7 rP (r\D)dr _ / °° rP(r, D)dr 

(1 + U)D - {r\D) - ^ p{rlD)dr ~ joo p(r , D)dr ■ (30) 

D is the 'raw' or derived distance to a source that is really at distance r. B is the bias in 
the derived distance D. The bias in the distance modulus will be /i^ as = 51og(l + B), and 
this is what is sought in this section. The one-sigma error bar for nuas is calculated with 
the usual higher moments of r\D. The expectation value of r given D is written as (r\D). 
P(r\D) is the conditional probability of r given D. P(r,D) is the joint probability which 
can be separated as 

P(r, D) = P(r)P(D\r) = r 2 n{r)P{D\r). (31) 

Here, P(r) is the probability of GRBs (observed or unobserved) being at a true distance of r, 
with this being r 2 n(r) where n(r) is the rate density of GRBs. The r 2 term accounts for the 
increased volume at greater distances. P(D\r) is the probability of measuring a distance D 
given that the burster is really at distance r. This probability contains all the information 
about the measurement uncertainties, the intrinsic scatter in the luminosity relations, the 
GRB luminosity function, and the gravitational lensing probabilities. 

The probability of measuring a distance D when the real distance is r, i.e., P(D\r), 
depends on many factors. In general it can be expressed as 



P(D\r) = [P M * (P L G)]P sample , 



(32) 
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for GRBs. Here the * symbol indicates a convolution. G is the Gaussian distribution of the 
measurement errors. Pl is from the burst luminosity function. The product PlG represents 
that there are more low-luminosity than high-luminosity bursts, so that deviations in that 
direction are more likely. Pm is the probability distribution of lensing with magnification M. 
The convolution in equation 32 produces the probability distribution of the burst brightness. 
The multiplication by P S am P ie is to account for the practicalities of detection of a GRB and 
its subsequent redshift measure as required for inclusion in this sample. 

The intrinsic scatter in the luminosity relations and the results of measurement errors 
can both be modeled together as a Gaussian distribution in log-space. This log-normal 
distribution is 

G(r, D) = e -^[(lnr-lnD)M nD ] 2 _ (33) 

The distance D is derived from the luminosity indicators (to get L) , the observed peak flux 
(Pbolo), and equation 23. The width of this distribution is given by 

a lnD = 0.5[K/L) 2 + (a Pbolo /P bolo ) 2 ] - 5 - (34) 

The luminosity L here is from the weighted average of the luminosities derived from the 
luminosity indicators (equations 12, 14, 16, 18, and 20). The uncertainty in that luminosity 
(ol) is from equations 13, 15, 17, 19, and 21 except that the systematic errors are set equal 
to zero. The reason is that the systematic errors likely arise from the lensing and Malmquist 
effects so it would count them twice to include them in the <jl values here. 

Firmani et al. (2004) derive the luminosity function for GRBs based on the peak 
flux distribution of BATSE bursts plus 33 GRBs with measured redshifts. The simplest 
acceptable models have the probability of a burst having a luminosity L as 

P L oc L- 1 - 57 . (35) 

The constant of proportionality does not matter as it divides out in equation 30. This can 
be translated (with equation 23) into a dependence on distance as 

P L (D)ocr~ 3A \ (36) 

for any individual burst. 

The lensing magnification M is the factor by which the unlensed brightness is multiplied 
to get the observed brightness of the burst. The average value of M must be unity. The 
majority of bursts will have M < 1 while there will be a small tail to high M. To take a 
specific example, for z — 3, the smallest value of M will be 0.81, the most likely value is 
0.89, the probability distribution has fallen to 10% of its peak at M — 0.82 and M = 1.16, 
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and there is a 1% chance of having M > 1.9 (Premadi & Martel 2005). The probability 
distribution of any particular magnification, P M , will depend on redshift. I have adopted 
the calculated P M curves out to z = 8 from Premadi & Martel (2005). 

The probability of including a GRB in my sample of 69 bursts (P S am P ie) depends on 
many factors which are poorly known and time variable. The detection of a particular burst 
requires some satellite experiment to be looking in the right direction at the right time to 
a threshold fainter than the burst (given its light curve shape). To a first approximation, 
the burst will be detected if its peak flux is brighter than the experimental threshold, which 
varies from satellite to satellite. Not only must the burst be detected, but it must also have a 
redshift determined from optical spectroscopy, and this depends on many variables relating 
to observing conditions and followup programs. To a first approximation, this probability 
depends on the brightness of the optical afterglow, which crudely scales as the peak flux 
of the burst. (Indeed, the Swift bursts have a probability of getting a redshift being only 
weakly dependent on the peak flux; with 28%, 38%, and 39% of the faintest-third, middle- 
third, and brightest-third, respectively, of the bursts discovered by Swift having redshifts.) 
So a reasonable approximation is to have P sa m P ie be only a function of Pboio- To quantify this, 
I have constructed a log N > p — log Pb i distribution for the GRBs in my sample, where N > p 
is the number of bursts brighter than any given Pboio- I have broken the sample into bursts 
according to the primary discovery satellite, as this will lump together GRBs discovered with 
a constant detection threshold and with similar dynamics in the optical followup community. 
With a knowledge of these thresholds and how followup efforts have changed over time, I have 
further grouped the BATSE and Konus discoveries together and well as the SAX, Integral, 
and HETE-2 events together. The resultant log iV>p — logP^ curves all have a steep 
region that is roughly -3/2 in slope which represents a uniform P sample above some threshold 
Pthreshoid- The curves all break at some peak flux and roll over in a manner well-described 
as a Gaussian. A Gaussian curve will have a long tail that must be cutoff somewhere, so I 
have adopted z = 7 as a cutoff due to the problems at getting spectra far enough into the 
infrared for higher redshifts. This probability function will be normalized to be unity for 
bright bursts, with this normalization choice dividing out in equation 30. Then, 



Here, W is a constant representing the sharpness of the cutoff in the discovery probabilities. 
For all satellites, I find that W = 0.5 fits well with the observed distributions. The Pthreshoid 
values I find are 10~ 5 erg cm" 2 s" 1 for BATSE and Konus bursts, 10" 6 erg cm -2 s" 1 for 
SAX, Integral, and HETE-2 bursts, and 3 x 10~ 7 erg cm" 2 s" 1 for Swift bursts. For use in 
equation 32, the peak flux has to be translated into a distance D for the given burst. As 




(37) 
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I am looking at brightness distributions for bursts in my sample, my formulation simulta- 
neously approximates the probabilities for both the burst detection as well as the redshift 
measurement. 

The number density (n(r)) of GRBs changes changes with distance, at least roughly 
proportional to the average high-mass star formation rate. Firmani et al. (2004) derive 
typical number densities as a function of redshift; 

n(z) oc e 2 - 3 * +c ^ 2 7(e 2 - 32 + 80), (38) 

where c = for z < 2 and c = 0.25 for z > 2. 

After all this, we can evaluate the bias with 

r 3 n(r)\P M * PLG]P S amviedr „ x 

(1 + B)D= J ° - . . L \ p sample . (39) 

J r 2 n{r)[P M * P L G\P sample dr 

This equation is to be applied to individual bursts. The required input is the satellite (to 
Set Pthre shoid), z (to select the Pm curve), L and Pb i (to determine D), as well as ctl and 
o Pholo (to determine a lnD ). 

Let me discuss the various competing effects to try to give some insight into the size 
and direction of the bias. There are two effects that will tend to make any given burst have 
a smaller measured distance than the GRB really is (i.e., D < r). The first is that there is 
greater volume at further distances so more GRBs will be included from distances greater 
than D than excluded from distances less than D. This is represented in equation 39 by the 
r 2 factor. The second is that GRBs are more common as we go to higher redshifts, so there 
are more GRB to be included with D > r than there are to be excluded with D < r. This 
is represented in equation 39 by the n(r) factor. There are also two effects that work in the 
opposite direction. The first is that if a burst appears farther than it really is (D > r), then 
it is less likely to be detected and a redshift measured. This is represented in equation 39 
by the P sa mpie factor. The second effect is that there are more GRBs with lower luminosity 
so that there will be more GRBs with D < r to get included in this sample than there are 
higher luminosity events to get excluded. This is represented in equation 39 by the P L factor. 
Thus we have two effects that are in competition with two other effects determining whether 
the bias will be positive or negative. The measurement error distribution (G in equation 39) 
has an average of zero, so this term does not introduce any bias. The lensing magnification 
distribution (P M in equation 39) has an mean of unity so that it also does not introduce any 
bias on average, even though it does provide scatter in the HD. 



The biases for my 69 burst sample have been calculated and presented in Table 7. 
Columns 1-3 are the usual GRB name, the detecting satellite, and the redshift. Column 4 
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is the peak luminosity (with units of erg s -1 ) as calculated from the luminosity indicators 
along with its uncertainty with no contribution from systematic errors. Column 5 is the 
bolometric peak flux and its uncertainty in units of erg cm -2 s _1 . Column 6 is the derived 
combined gravitational lensing and Malmquist biases expressed as a shift in the distance 
modulus, fiuas and its uncertainty in units of magnitudes. 

In the absence of gravitational lensing, accurately measured distances (i.e., with small 
ctl and &p bolo ) will result in very small biases. Indeed, the bias will scale as the square of 
°"inD (e.g., Gonzalez & Faber 1997). For GRBs, a^D is not small, so the first expectation 
is that the biases will be quite large. But a glance down the last column in Table 7 shows 
that the biases are quite small. The average value of fiuas is 0.03 mag while it has an RMS 
scatter of 0.14 mag. Thus the biases are small compared to the the typical error bars in /i 
as well as the typical changes in \i for the model HD at moderate redshifts. There are two 
reasons why the GRB biases are small. 

The first reason for small biases is that the bursts with redshifts are on average the 
same ratio above the experimental threshold. That is, the faint bursts were discovered with 
satellites with low thresholds while the bright bursts were discovered with experiments with 
high thresholds. Thus the median values of Pboio/ Pthreshoid is 0.87 for BATSE and Konus, 
0.52 for SAX, HETE-2, and Integral, and 1.08 for Swift - despite the thresholds varying by 
over a factor of 30. This results in the values of Pboio/ Pthreshoid having no significant trend 
with redshift; with median values of 1.09 from < z < 1, 1.27 from 1 < z < 3, and 0.64 
from 3 < z < 7 - despite the Pboio changing by a factor of 1200. That Pboio/ Pthreshoid is 
largely independent of redshift causes the bias corrections to also be largely independent of 
redshift, and thus that the shape of the HD does not change. 

This first reason for the low biases for GRBs is identical to one of the reasons why 
the Malmquist bias is claimed to be negligible for supernovae in the HD (Perlmutter et al. 
1997; 1999; Knop et al. 2003). In this case, the low redshift sample of supernovae (from 
Hamuy et al. 1996) had relatively bright peak magnitudes yet the discoveries were made 
with a relatively high detection threshold, while the high redshift sample supernovae had 
relatively faint peak magnitudes yet the discoveries were made with a relatively low detection 
threshold. The logic was that both samples had similar brightnesses above thresholds at the 
time of discovery so that any Malmquist bias would shift the two samples by similar amounts 
resulting in no change of shape to the HD. 

The second reason for the small biases for GRBs is that the product r 2 n(r)P L P samp i e 
is largely flat for most of the GRBs in my sample. That is, fortuitously, most bursts have 
that product roughly constant across the likely range of distances. In this case, the integrals 
in equation 39 will simply return D, so that both B and jibias will be near zero. Let me 
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illustrate this with a specific case of how the product changes from z = 1 to z = 2. The 
r 2 term increases by a factor of 5.5, the n(r) term increases by a factor of 16, the Pl term 
decreases by a factor 0.24, and the P sa mpie term decreases by a factor of 0.21 for GRBs farther 
than the threshold. In all, the product r 2 n(r)PLP samp i e changes by a factor of 0.44 across 
this range in redshift. And most bursts have their likely range of redshift being much smaller 
than this, so that the change in the product across their likely distances will generally be 
much smaller resulting in a small Huas- 

This second reason is similar to the case for small Malmquist effects for the supernova 
HD (Perlmutter et al. 1997; 1999; Knop et al. 2003). These papers explicitly cite the 
offsetting effects that cancel out as part of their justification that the Malmquist bias can be 
neglected. 

In all, we see that the gravitational lensing and Malmquist biases are much smaller 
than the quoted error bars and smaller than typical variations in the models. Even though 
the systematic biases are low, there is still substantial scatter introduced into the GRB HD 
by these two effects. At a redshift of 3, the RMS scatter in /j, due to variations in the 
lensing magnification is roughly a quarter of a magnitude (Premadi & Martel 2005), and 
this is comparable to the systematic scatter in the luminosity relations. Thus, I attribute 
the 'extra' scatter to the lensing and Malmquist effects. 



To put constraints on parameters in cosmological models, we must first adopt a set of 
priors as constraints on allowed values for parameters and we must identify which parameters 
are statistical 'nuisance parameters' (whose value we do not care about) for marginalization. 

Different workers will choose different priors. Throughout this paper, I will take the 
universe to be flat (Qm + = 1). For the cases where I am examining possible variations 
in the Dark Energy, I will adopt priors of either < Vt M < 1 or 0.23 < Q M < 0.31. 

The identity of nuisance parameters depends on the science question being asked. If 
we are interested in constraining the cosmological parameters, then we really do not care 
about the value for the slope of the V — L relation (or any of the other nine fitted constants 
in the luminosity relations). The process for marginalizing these nuisance parameters is to 
integrate over them. For example, if we assumed a flat universe and w — — 1 and are only 
interested in Qm, then the probability for a particular value of VLm is given by 



7. Cosmological Constraints 




(40) 
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where a\ and b\ are the intercept and slope of the first luminosity relation, and so on. These 
probabilities are related to the observed chi-squares as 

P oc e-* 2 / 2 . (41) 

If we are interested only in the Dark Energy parameters, then we might consider VLm as a 
nuisance parameter that also needs to be marginalized. 

What are the best fit parameters for the cosmological models? Here, I will examine 
three separate questions: 

If Dark Energy is constant with w — — 1 and the universe is flat, then what is the 
best value for Vl M from GRBs alone? This is somewhat of a test case, as supernovae and 
other evidence already tell us that Qm = 0.27 ± 0.04 in this case. Also, other GRB studies 
have results that can be used for this same question. Thus, Firmani et al. (2006) find 
Q M = 0.291^ while Xu, Dai, & Liang (2005) find Q M = 0.30±g;g|. Figure 10 presents 
P(VIm) with an arbitrary normalization. I find that VLm = 0.39^g'og. 

What constraints can we place on the Dark Energy parameters (wq and w a ) if we are 
not interested in the value of Qm in a flat universe? For this question, we have to marginalize 
over the luminosity relation parameters as well as Qm over the range 0-1. Figure 11 shows 
the resulting one-sigma and two-sigma confidence regions. The best fit set of parameters 
is wo = 0.1 and w a = —0.9. If we are only interested in one of these parameters, then the 
one-sigma uncertainty range, then this will be spanned by the one-sigma region in Figure 
11. Thus, w = 0.1^5 and w a = — 0.9^^2 at the one-sigma level and w = O.ll^J and 
w a = — 0-9-f\ at the two-sigma level. But these simple one dimensional error bars do not 
represent the correlation in Figure 11. So perhaps a better representation of the one-sigma 
confidence region would be 

w a + 2.2w = -0.68l°;l, w a - 0A5w = -0.95l^. (42) 

In the thin dimension, the width of the one-sigma region is 0.73. 

In this case, the Cosmological Constant (w = — 1) has a marginalized probability that 
is only 12% of the peak probability, and this puts that position just outside the two-sigma 
confidence region. However, the best fit area on this diagram corresponds to a relatively 
shallow slope for z < 0.5, and this is inconsistent with the supernovae data. If the supernova 
data is used (either as a prior or as part of a joint fit) to exclude much of the region sub- 
stantially below and to the right of the Cosmological Constant point, then the Cosmological 
Constant point will likely go inside the one-sigma confidence contour. Another reason to 
not reject the Cosmological Constant from this analysis is that the chi-square for the best 
fit luminosity relations (i.e. without any marginalization) is 70.7 (for 66 degrees of freedom) 
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for the best point in Figure 12 and is 72.3 (for 69 degrees of freedom) for the Cosmologi- 
cal Constant point. Naturally, the use of three additional fit parameters has improved the 
chi-square, but an F-Teast shows that the improvement is much too small to justify the use 
of the additional parameters. Thus, the improvement in chi-square is too small to reject 
the Cosmo logical Constant. Another way to see that the Cosmological Constant should not 
be rejected is simply to recall that the reduced chi-square for its HD fit is 1.05, and this 
indicates an acceptable model. 

If we accept the prior work yielding Qm = 0.27 ± 0.04 in a flat universe, then what are 
the limits on the Dark Energy parameters? For this question, we have to marginalize over 
the luminosity relation parameters as well as Qm over the range 0.24 - 0.31. Figure 12 shows 
the resulting one-sigma and two-sigma confidence regions. The case with a closely restricted 
Vl M has a best fit location of w = 0.2±0.6 and w a = — 1.4 ± 1.2 (at the one-sigma confidence 
level), while its unmarginalized chi-square is 70.7. The one-sigma confidence region can be 
represented by 

w a + 2.1w = -1.0 ±0.3, w a - 0.48wo = -1.5 ±1.5. (43) 

In the thin dimension, the width of the one-sigma region is 0.6. The comparison between 
Figures 11 and 12 shows that they are largely the same. Nevertheless, the two-sigma region 
for Figure 11 is somewhat larger than in Figure 12, as appropriate for its looser constraint 
on Q M . 

In the case with the more restrictive prior on M , the Cosmological Constant (w = — 1) 
has a marginalized probability that is only 11% of the peak probability, and this puts that 
position just outside the two-sigma confidence region. Again, the results of this fit do not 
allow for the rejection of the Cosmological Constant. 

In both Figures 11 and 12, we find the Cosmological Constant point being just outside 
the two-sigma confidence region. But we cannot reject the Cosmological Constant on this 
basis, as discussed above. The real bottom line is that the simple concordance case (w — — 1 
and Qm — 0.27 in a flat universe) provides a nice fit to the data (see Figure 7) as evidenced 
by a fine reduced chi-square of 1.05. By allowing three additional parameters to vary, we 
can always get a slightly better reduced chi-square, in this case with an improvement from 
72.3 to 70.7. In all, I find that the GRB HD is fully consistent with a constant Cosmological 
Constant. 
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8. How Robust is the GRB HD? 

How robust are the results coming from the fits to the GRB HD? That is, for reasonable 
differences in various choices in the analysis, will the chi-square contours (as in Figures 11 
and 12) change significantly? I can perform a variety of 'sanity checks' and I can make a 
variety of procedural changes to see the sensitivity of the results to the changes. In this 
section, I will try seeking for systematic changes in the HD and fits as various selections are 
varied. 

The results are tabulated in Table 8, with the various cases specified and the number of 
GRBs used in each case. To quantify the confidence contours, I can determine the minimum 
chi-square, the position of the long axis and the width perpendicular to that axis (expressed 
as a one-sigma error bar on that position), the acceptable position along that axis (expressed 
as a range of wq), the chi-square for the concordance model, and the chi-square for the 'Riess' 
model. The first line in Table 8 is for the 'standard' data set from Tables 2-4, corresponding 
to the case displayed in Figure 12. Subsequent lines are for different cases and should be 
compared back to the first line to see whether the differences are significant. 

The HD has only two GRBs with z > 5, so we can wonder how sensitive the results are to 
these two points? That is, are the conclusions dominated by just two bursts, or did biases in 
these two events significantly shift the cosmological constraints? To address these questions, 
I have repeated the entire fitting of the GRB HD (as in Figure 12) while deleting GRB 
050904 and GRB 060116 from the data set. We see the that the chi-squares are somewhat 
smaller due to eliminating two bursts (as expected). Also, the chi-square differences between 
the various models is lower, as expected when the two points with the highest redshifts are 
arbitrarily ignored. The position of the confidence contours are essentially unchanged. In 
all, the constraints on cosmology and the conclusions are not affected by a deletion of the 
two highest-redshift bursts. 

How sensitive are my cosmological fits to the rejection of the outliers as itemized in 
Section 3. Recall, that I rejected three bursts and two rise times as being outliers. One way 
to test this is to repeat the fitting to the GRB HD except to include all the rejected input. 
The properties of the new chi-square contours are tabulated in Table 8. We see that the 
minimum chi-square is very large (209.9 versus 70.7 with the outliers rejected), but this is 
exactly as expected. Most of this large increase comes from just the one burst GRB 980425. 
The central axis of the long-thin confidence contours is the same, although the high position 
on the HD and low redshift of GRB 980425 pushes the minimum chi-square to w > 2 along 
this axis. This is not surprising, as any one very distant outlier can force a substantial tilt 
in the best fit and hence push the contours around. If we include the data for four of the 
outliers not including GRB 980425 (i.e., GRBs 020819 and 050315 as well as the rise times 
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for GRBs 990123 and 030328), then we get the contour as shown in Table 8. The chi-squares 
are still high when compared to the standard case (as expected). But largely we see that 
the position of the contours and the relative chi-squares are not affected by the inclusion of 
the four outliers. 

In Section 2, I rejected five GRBs which have reported "redshifts of too low a confidence 
to be used", so we can ask whether these rejections are reasonable? For GRB 980326, the 
reported redshift of z ~ 1 (Bloom et al. 1999) was solely based on interpreting a marginally- 
detected bump in the light curve as being a SN1998bw-like supernova. However, we now 
know that the GRB/supernova events span a wide range of peak absolute magnitudes, so 
the derived redshift must have a very large uncertainty, even assuming that the bump arises 
from a supernova. And the z ~ 1 claim is dubious since the host galaxy would then have 
to be at least seven magnitudes fainter than L* and this is very unlikely (Fruchter et al. 
2001a). So the reported redshift is much too low in confidence for GRB 980326 to be used 
in a Hubble diagram. GRB 011121 has a reported redshift of z — 0.362 on the basis of a 
nearby galaxy with no connection to the optical transient position (Garnavich et al. 2003). 
Rather, with EST imaging, the host appears to be a separate 'blue knot' of unknown redshift 
(Bloom & Price 2002). GRB 020305 has a galaxy on top of the optical transient position, 
for which a low signal-to-noise spectrum yields proposed redshifts of z ~< 2.8, z ~ 0.2, 
z > 0.5, and z = 2.5 (Gorosabel et al. 2005). Clearly, any of these four possibilities has 
too low a confidence for GRB 020305 to be used in any Hubble diagram. GRB 050803 has 
only an x-ray error circle, and this has "several faint sources in the uncertainty circle", for 
which the brightest happens to have a redshift of 0.422 (Bloom et al. 2005). An error circle 
with 12" diameter will always have some brightest source and any connection to the GRB 
can only be conjectural. GRB 060123 also only has an x-ray error circle and naturally there 
is some brightest galaxy in that circle (Butler & Bloom 2006). A long spectrum from the 
Gemini 8-m telescope shows only one certain emission line, for which redshifts of z = 1.099, 
z = 0.193, and z = 0.562 have been proposed (Berger et al. 2006), so that all we have is 
a galaxy of unknown redshift and unknown connection to the GRB. In all, the five rejected 
GRBs have suggested redshifts that are of much too low a confidence for inclusion in any 
HD work. 

A related question is whether any of my 69 GRBs have a redshift of sufficiently low 
confidence that some would choose to reject those GRBs? This is partly a difficult ques- 
tion because the available information is very scattered in nature and hard to lead to a 
quantitative confidence level. Thus, we cannot simply perform fits as a function of some 
redshift-confidence threshold. Instead, I can perform the cosmology fits where I eliminate 
any of my 69 GRBs for which the redshift has been seriously questioned. I find that only 
three of my GRBs have had serious questions raised about their quoted redshifts; GRBs 
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050802, 060108, and 060116. GRB 050802 has a spectrum of its optical afterglow with sev- 
eral absorption lines indicating a redshift of 1.71 (Fynbo et al. 2005c), however, McGowan et 
al. (2005) place a limit of z < 1.2 based on a UVOT detection in the far ultraviolet. But the 
optical redshift is of very high confidence with a good spectrum and many lines (J. Fynbo 
2005 private communication), while the redshift of z = 1.71 is actually consistent with the 
UVOT detection (K. McGowan 2005, private communication). As such, I have confidently 
included this burst into my HD. GRB 060108 has its redshift based on a fit to the spectral 
energy distribution of the optical afterglow, with the best value being z = 2.03 and an upper 
limit being quoted as z < 2.7 (Malendri et al. 2006). As such, perhaps the uncertainty 
in the redshift is too large for inclusion in the GRB HD? GRB 060116 has a well-observed 
optical/infrared afterglow whose spectral energy distribution yields z = 6.60±0.15 (Grazian 
et al. 2006). A more detailed analysis of the same data (Piranomonte et al. 2006) returns 
the same redshift, but another possibility of 3.8 < z < 4.5 is presented with a lower likeli- 
hood. The two possibilities are distinguished by either low extinction (for z = 6.60) or high 
extinction (for z ~ 4). Piranomonte et al. (2006) report a "very poor" late spectrum with 
possible light down to 700 nm wavelength, but this is contradicted by the lack of a detection 
with the HST ACS/F775W filter (Tanvir et al. 2006), so that the z = 6.60 value is strongly 
preferred. Despite having reasonable confidence for the inclusion of GRBs 050802, 060108, 
and 060116, we can ask what would be the effect on the cosmology fits if these three bursts 
are excluded from the analysis. The resultant chi-square contours are presented in Table 8. 
We see that the rejection of the three 'questionable' redshift bursts does not significantly 
change the cosmology fit results. In all, the GRB HD appears to be robust to changes in the 
confidence level for the measured redshifts. 

In principle, the weighted average of the distance moduli from all available luminosity 
relations will combine all information and produce the best value. But we can worry whether 
one or more of the relations is adding in more noise than signal. To test this, I have performed 
all the fits again while successively eliminating each of the luminosity relations. The chi- 
square contour information is in Table 8. As expected, the minimum chi-square varies around 
a bit. However, the primary axis of the confidence region and the relative chi-square values 
vary little. An examination of the HDs reveals a well formed set of points following along 
a smooth curve with no signs of systematic shifts. As such, I see no evidence that any one 
relation is adding more noise than signal. In another trial, I have constructed a HD (see 
Figure 13) from the standard data set with the elimination of the distance moduli derived 
from the two relations (V — L and trt — L) with the most scatter in the calibration curves. 
The plot points fit smoothly along the model curve and all looks fine. Figure 13 can be 
directly compared with Figure 7. As expected with the two relations carrying 32% of the 
total statistical weight, the typical error bars in Figure 13 are substantially larger than those 
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in Figure 7. In all, there is no evidence that any of the relations is doing more harm than 
good, so it is best to use all the relations and take advantage of all the available distance 
information. 

If each relation is taken one at a time, then are the resulting chi-square contours consis- 
tent with the contours from all five relations? To test this, I have constructed HDs for each 
luminosity relation alone. As expected, the scatter is much larger than displayed in Figure 7. 
The contours are presented in Table 8. The contours are consistent with the standard case, 
with the expected shifts as will always be seen when looking at a subset of the data. The 
largest shift is for considering the T\ ag — L relation. The E pea k — E 1 relation is surprisingly 
broad, although it is still consistent with the five-relation contours at the one-sigma level. 
In all, I see nothing that points to any one relation being inconsistent with the final result. 

One possible selection effect could arise from the diversity of GRB spectral shapes, 
especially at higher energies. A high redshift population of GRBs will have their higher 
photon energies shifted to the detector's bandpass and then be selected for those with harder 
spectra (i.e., with less negative (3 values). In this way, the detected bursts at higher redshifts 
might have a harder spectrum than the sample of detected bursts at low redshift, even in a 
case with no evolution of the full GRB population. I have three reasons for thinking that 
such effects are minimal: (A) The high redshift bursts are necessarily of high luminosity and 
hence (by the E pea k — L relation) have a high intrinsic E pea k. Schaefer (2003b) has shown 
that the redshifting comes close to exactly offsetting the intrinsic E pea ^ so that bursts at 
all redshifts will have the same apparent E pea ^. (This is actually a function of the observed 
peak flux, and the insensitivity of of the measured E pea k on redshift is the reason that the 
observed E pea k distribution is allowed to be as narrow as it is.) In this case, the effects of 
diversity in the GRB spectra on the burst triggering will be identical at both high and low 
redshifts. (B) The high energy band contributes only a relatively small fraction of the flux 
in the trigger range, so variations of (3 will have little effect on burst detection. To take a 
specific example, for the usual BATSE trigger range of 50-300 keV, a burst with the median 
E pe ak = 300 keV (Mallozzi et al. 1995) will have zero effects on the trigger caused by varying 
f3 since the spectral break is above the trigger band. If E pea k = 100 keV, then a variation of 
—2 < (3 < —3 will only cause the counts in the BATSE trigger range to change by 6%. For 
very low values of E peak , say at 30 keV, the effect of varying (3 over such a large range can 
produce a factor of two variation in the counts in the trigger band. For the Swift detectors, 
the trigger band is 15-150 keV, hence making the effects of varying (3 zero for most bursts 
and small for most of the remainder. In all, the triggering effect of (3 variations is minimal. 
(C) If a selection effect on (3 as a function of redshift exists, then we might expect to see a 
correlation between f3 and z for the detected bursts. But such is not seen from bursts with 
measured (3 values in Table 2. For example, the weighted averages of (3 for z < 1 and z > 2 
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is -2.31 ± 0.05 and -2.26 ± 0.06 respectively. 

A means of pointing to the improvement implemented by the luminosity indicators is to 
construct a HD with the bursts artificially set to the same luminosity. This is displayed in 
Figure 14 for the concordance cosmology. This figure shows a horrible match with immense 
scatter. As such, the comparison between Figures 7 and 14 demonstrates that the luminosity 
relations are contributing vital distance information. 

In summary the GRB HD and the fits to cosmology appear to be robust. 

9. Comparison with Supernovae 

For Hubble diagram work, a comparison of GRBs and supernovae is inevitable. 

The great advantage of supernovae for the HD is that they are substantially more 
accurate standard candles than are GRBs. This will translate into better constraints on 
cosmological parameters. But this advantage for supernovae is not as large as some might 
expect. The one-sigma scatter of supernovae about the concordance cosmology in the HD 
is 0.22 mag for the best 31 supernovae out of 42 or 0.36 mag for all 42 (Perlmutter et al. 
1999), 0.27 mag for 0.1 < z < 1 and 0.29 for 1 < z even after the selection of the very best 
supernovae (Riess et al. 2004), and >0.25 mag (Astier et al. 2005). It is clear that the real 
accuracy of supernovae is ~ 0.25 mag only after much selection, so the real accuracy without 
a posteriori selection (other than throwing out > 3 — o outliers) is more like 0.3 mag or 0.35 
mag. This is to be compared to the residuals of my GRB HD where the RMS scatter is 0.65 
mag. So the comparison between supernovae and GRBs is that the scatter in the residuals 
around the concordance cosmology is 0.3-0.35 mag to 0.65 mag. That is, supernovae are 
roughly a factor of two more accurate than GRBs for HD work. This factor of two is likely 
to be substantially smaller than many people would expect. And there is every prospect 
of significantly improving the GRB luminosity indicators (see Section 10), so this factor is 
likely to get substantially smaller over the next few years. 

A great advantage of GRBs for HD work is that they go out much further in redshift 
than the many ground-based supernovae (z < 1) or space-based discoveries from either HST 
or SNAP (z < 1.7). GRBs are already numerous from 1.7 < z < 6.6, and it might be 
possible in the future to get an occasional GRB out to even z = 10 (Lamb & Reichart 2000; 
Fenimore & Ramirez- Ruiz 2000; Schaefer, Deng, & Band 2001; Bromm & Loeb 2002). The 
use of high redshift GRBs will provide a long lever-arm for measuring changes in the slope 
of the HD. That is, model HDs will have small differences in predicted distance moduli for 
z ~ 1 yet will have substantially larger differences at high redshift (see Figure 15). For 
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example, the differences between the concordance cosmology and Riess' best fit cosmology 
is -0.03 mag at z — 0.5, 0.04 mag at z = 1, 0.44 mag at z — 3, and 0.83 mag at z — 5. More 
to the point for a comparison between GRBs and supernovae, the difference for the highest 
redshift supernova (at z = 1.7) is only 0.15 mag while the difference for the current highest 
redshift GRB (at z = 6.6) is 1.00 mag. Thus GRBs have a large advantage by using the long 
lever arm of high redshifts since the models predict large differences at high redshifts. 

A common argument is that the HD is uninteresting for redshifts higher than the SNAP 
range because the universe will be matter dominated and so the Dark Energy will be too 
small to affect the universal expansion at high redshifts. But this argument is wrong because 
it ignores the long lever arm of going to high redshifts in yielding the slope of the HD. To 
take a specific example, suppose that we wanted to measure the slope of the HD above z — 1. 
Would it be better to measure one supernova (with er M = 0.30 mag) at z — 1.7 or one GRB 
(with (j n = 0.65) at z — 6.6? The model differences between the concordance and 'Riess' 
cosmologies is 0.15 mag at z — 1.7 so the accuracy from a single supernova is twice the size 
of the effect being sought. Alternatively, the model differences at z — 6.6 is 1.00 mag, so 
the single high redshift GRB will detect the differences at close to the 1.5-sigma level. In 
this case, the single GRB is three times better than the single maximal-redshift supernova. 
That is, the long lever arm for GRBs more than makes up for their lower accuracy. Or 
another way of looking at it is that a single GRB will provide more information than a single 
supernova. 

The above common argument is also wrong at a higher level. This argument declares 
the 1.7 < z < 6.6 region to be uninteresting within the concordance cosmology. But this is 
presumptive that the concordance cosmology is correct. (Similarly, back in 1997, I was told 
by a theorist that it was pointless to measure the HD to z — 1 since it would merely confirm 
the then concordance cosmology of VL M = 1 and = 0.) Science advances by exploring 
unexplored regions and by performing critical tests of standard wisdom. (And the standard 
wisdom of the now concordance cosmology is only a few years old.) Who knows what we 
will find in the 1.7 < z < 6.6 HD, and we won't know unless we look. It would be unwise 
for the community to ignore 1.7 < z < 6.6 as uninteresting, especially as the GRB data is 
currently flowing in for free from Swift. 

A disadvantage of supernovae for HD work is that the optical light can be dimmed 
with distance. Two serious mechanisms have been proposed. The first is that the optical 
light will be achromatically dimmed by grey dust and hence not corrected for in the usual 
dereddening (Aguirre 1999a; 1999b). Simple versions of this mechanism have been excluded, 
but the possibility that the dust density changes with redshift has not been excluded (Riess 
et al. 2004). The possibility of grey dust changing over time to match the predictions of 
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the concordance cosmology seems contrived, but it is possible. The second mechanism is 
that the optical light will be achromatically dimmed by refraction in Lyman-a clouds along 
the line of sight (Schild & Dekker 2005). This idea remains largely unexplored. GRBs as 
cosmological tools are completely immune to both problems. 

A likely serious problem for the supernova HD is that the progenitor population of 
Type la supernova undoubtedly evolves in time, so that high redshift events might have a 
substantially different calibration from nearby events, thus leading to a distortion of the HD. 
This is a reasonable possibility because the progenitors formed in the young universe will 
have lower metallicity than progenitors formed in recent times, and the metallicity might 
have a noticeable effect on the supernova peak brightnesses and decline rates. The primary 
defense is that the spectra of supernovae at z < 1 appear similar to those of nearby events 
(Perlmutter et al. 1997; 1999), and this is good for putting some crude upper limit on the 
size of the evolution effects. However, the peak brightnesses of supernovae are correlated 
with the galaxy type, and galaxy types evolve as we look back to higher redshifts, so we have 
observational evidence that the evolution effect is significant. Detailed calculations under 
various scenarios (Dommguez, Hoflich, & Straniero 2001; 2003) show that the progenitor 
evolution effects are roughly 0.2 mag from z — to z — 1. This is comparable to the size 
of the cosmological effects over the same redshift range. As such, the supernova HD cannot 
be used for precision cosmology until the evolution issue is resolved. And there is currently 
no resolution of the long-standing question of whether the progenitors are recurrent novae, 
double white dwarfs, supersoft binary systems, or symbiotic stars. Without knowing the 
identity of the progenitors, how can we evaluate the evolution of peak brightnesses as we 
go to high redshift? In all, supernova cannot be applied to HD questions where systematic 
changes of less than ~0.2 mag from low to high redshift are critical. 

Do GRBs have the same problem with evolution? I argue that GRBs do not have 
any problems with evolution. My reason is that the physical mechanisms that create the 
luminosity relations are thought to be based on light travel times, the degree of relativistic 
beaming, and energy conservation in the shocked material - with none of these mechanisms 
changing as we look back in our universe. The metallicity of the GRB progenitor or the 
surrounding ISM does not change the speed of light or the relativistic effects that are the 
basis for the luminosity relations. The metallicity might affect the GRB luminosity, but then 
the physics of the luminosity relations would simply indicate the corresponding luminosity 
and the distance would be correctly deduced. The population of GRBs might drift, for 
example, in luminosity as we look back to high redshift, but the luminosity indicators will 
still give the correct luminosity for each individual burst. And this is all we need for GRBs to 
have zero evolution effects. This argument should be examined further, but in the meantime, 
we are left with a situation where GRBs have arguably zero evolution effects whilst supernova 
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have currently unknown evolution effects that can get as large as the signal being sought. 

Do either GRBs or Type la supernovae have an advantage concerning the number of 
outliers that must be rejected? In Section 3, I rejected three GRBs and two rise times out 
of 72 bursts. For comparison, Perlmutter et al. (1999) presented their main result with 6 
rejected outliers out of 60 supernovae. Riess et al. (2004) rejected 29 supernovae (out of 186) 
to form their 'gold sample', with the causes largely being due to uncertain classifications, 
too few points in the light curves, and high extinctions. Many well-observed supernovae 
are known to be distant outliers (not even counting the super-luminous SN 1991T and sub- 
luminous SN 1991bg events); including S Andromedae (de Vaucouleurs & Corwin 1985; 
Fesen, Hamilton, & Saken 1989), SNLS-03D3bb (Howell et al. 2006), SN 2005hk (Phillips et 
al. 2006; Stanishev et al. 2006), SN 2002cx (Li et al. 2003), SN 2003gq (Jha et al. 2006), SN 
2005P (Jha et al. 2006), SN 2005cc (Antilogus et al. 2005), and SN 1006 (Schaefer 1996a; 
1996b). Apparently, Type la supernovae have a higher percentage of outliers than GRBs, 
but I do not think that this is important for the relative merits. Outliers are not important 
for either Type la supernovae or GRB HD work because the outliers can be readily identified 
and rejected for both types of explosions and because the outliers occupy a suitably small 
fraction of the overall event population. 

Supernovae have a substantial advantage for HD work because their physics is well 
known, whereas the physics of GRBs is relatively much less well known. That is, the mech- 
anisms of the supernova explosion have been very closely studied by many groups and we 
have excellent and extensive observational material. In contrast, the basic scenario for GRBs 
(internal shocks within a relativistic jet produced by a core collapse of a very massive ro- 
tating star) is only a few years old, and many of the details are still being debated. This 
situation leaves a comforting feeling when dealing with the supernova HD; while we are left 
to ponder whether there are any now-hidden surprises in the GRB HD. This advantage for 
supernovae might only be psychological, however, because the calibrations and corrections 
in all supernova cosmology projects are entirely empirical. That is, there is zero contribution 
from supernova theory in deriving the HD. So nowhere has the advantage been used. The 
GRB HD is also entirely empirical (see Figures 1-5), and thus it is operating on the same 
basis as the supernova HD. 

I have just made a strong case that the GRB HD has some good advantages over the 
supernova HD and that its disadvantages are not as bad as some might expect. So which is 
'better', or more to the point, which should get attention and resources? For this, the answer 
is obvious that both should get attention and resources. The reason is simply that no one 
method can convince the community by itself, so that multiple methods giving concordant 
results are required. (This was a lesson learned deeply with the original supernova HD work 
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in the late 1990s.) That is, any answer from the GRB HD alone will not convince the 
community, just as any answer from SNAP alone will not convince the community. What is 
needed is completely independent methods that concur. To this end, our community must 
strongly support the SNAP mission, just as GRB HD efforts should be supported. 

10. Future Work 

An obvious extension of this work is to perform simultaneous fits with GRBs and su- 
pernovae. Care would have to be taken that each data set is handled correctly, for example 
with the GRB luminosity relations being fitted for every set of trial cosmological parameters. 
The supernovae would dominate the delineation of the HD for 1 < z while the GRBs would 
provide all the information for the HD for 1.7 < z < 6.6. This use of all data will provide 
the best possible constraints on the shape of the Hubble diagram. 

Another obvious extension of this work is to test a more versatile expansion for w(z). 
The w(z) = Wq + w'z expansion has problems with divergences for z > 1. The w(z) = 
Wo + w a z(l + z)~ l expansion makes a presumption that most of the change in w occurs at 
z ~ 1 (which was about the only place where it could be tested with the supernovae) which 
might or might not be appropriate. Perhaps an expansion like w(z) = w + w a z[z t + z)^ 1 , 
where z t is the redshift of the transition (Rappetti, Allen, & Weller. 2005). With the 
supernova plus GRB data covering such a wide range of redshift, it might be possible to go 
to additional parameters. 

I can think of a variety of places in which substantial improvements can be made within 
the next year or two. The most obvious of these is that the luminosity relations can be 
improved. To this end, Firmani et al. (2005; 2006) have found an improvement in the E peak — 
L relation by including duration information, and Li & Paczyhski (2006) have improved the 
Vl relation with an alternative normalization factor and a different method of smoothing 
the light curve. Other improvements can be expected, and let me list here a few ideas: 
(1) We redshift and dilate the various luminosity indicators to the rest frame of the GRB 
with factors of 1 + z, but we measure the luminosity indicators and peak fluxes from light 
curves created for energy bands fixed in the Earth rest frame. With Swift data and analysis 
software being public, it is now possible for investigators to calculate the luminosity indicators 
with light curves for standard energy bands in the rest frames of the individual GRBs. (2) 
The proliferating definitions of variability all have many arbitrary choices, so I wonder if a 
completely different approach might result in a tighter relation. One idea would be to look 
at the slope or intercept of the Fourier power spectrum, but I expect that this will be poor 
since the Fourier spectra are more appropriate for long data streams. A better idea is to use 
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a wavelet activity spectrum (as in Walker, Schaefer, & Fenimore 2000) as wavelets are local 
measures and the activity directly measures the amplitude and time scale of the variations. 
(3) Firmani et al. (2005; 2006) use a parameter called T 0A5 , which is the duration during 
which the burst light curve is brighter than 45% of the peak flux. The choice of the 45% 
level is only left over as a historical artifact and there is no real likelihood that this level 
is optimal. An obvious study is to vary this level and seek the minimal scatter about the 
calibration equation. And other duration definitions might be explored, for example T 90 , 
Tgo/Apeafc, T 0A5 /N peak , or even trt- (4) We now have many measurable parameters that 
are correlated with luminosity, and these might merely be projections of higher dimensional 
correlations. It might be profitable to perform a multivariate analysis to try and discover 
some sort of a 'fundamental plane' which can be used to substantially reduce the scatter in 
the luminosity relations. 

Another substantial improvement in the near future will be that Swift (along with 
ground based observers) will be providing more bursts to include in the GRB HD. The 
Konus-Wind and Suzaku-WAM experiments will provide tight limits on the E pea k values for 
many bursts. Swift has produced 22 GRBs for my sample in the year before the last entry. 
In two years, Swift will have a total sample size of roughly 75 bursts. One advantage of such 
a large Swift sample is that they will all come from a single instrument. Such a homogeneity 
will produce more consistent luminosity relations and ease the calculations of Malmquist 
bias. The sample over the next few years will also be important as an independent data 
set from the one in Table 2, and this will be a conceptually important means of testing the 
results from the current sample. 

It is possible to design a specialized space mission that will greatly improve the accuracy 
of the GRB HD. For example, Lamb et al. (2005) describe a possible MIDEX-class mission 
that will measure z, E peak , S, and T jet for >800 GRBs with 0.1 < z < 10 during a two-year 
mission. Such a program would measure VL M and w to an accuracy of roughly ±0.06. 

The future will have the details of the GRB HD examined closely by many groups, with 
all the normal problems being raised and resolved. With the GRB data now pouring in 
for free from Swift, our community will be seeing a lot of HD work over the near future. 
But ultimately, any conclusions from the GRB HD will have to be tested and confirmed by 
independent methods. 
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11. Conclusions 

In this paper, I have constructed a Hubble diagram for 69 Gamma-Ray Bursts with 
redshifts from 0.17 to >6. The basis for evaluating the luminosity distances is five luminosity 
indicators (r/ aff , V, E peak , tj et , and r RT ). The construction of the GRB HD has simply followed 
in the footsteps of many prior papers, and no new techniques are being used here. What is 
new is that I am using over one order of magnitude more data than any previous work. This 
is achieved by the simultaneous use of all five indicators (rather than only one or two at a 
time) as well as the use of recent bursts discovered by Swift. My conclusions are: 

(1) The GRB Hubble diagram defines a well-behaved curve with half the bursts at 
z > 1.7 and seven GRBs with z > 4. This bodes well on the utility and accuracy of 
cosmological information from the GRB HD. 

(2) The average one-sigma uncertainty for a single GRB is 0.65 mag in the distance 
modulus. This is a factor of two larger than for the average one-sigma uncertainty for a 
single supernova. This bodes well for the future, as there are a variety of means by which 
the GRB luminosity relations can likely be substantially improved. 

(3) The gravitational lensing and Malmquist biases are small, with the average bias of 
0.03 mag and the RMS scatter of this bias being 0.14 mag. This surprisingly small bias 
arises from both of two reasons; first that the average burst brightness (compared to the 
experimental threshold) varies by less than a factor of two as redshift goes from low to high, 
and second that the various competing effects (from increasing volume, increasing number 
density, decreasing luminosity function, and decreasing detection probability with increasing 
distances) all come close to canceling each other out for typical burst properties. 

(4) The 1.7 < z < 6.6 region of the HD is of high interest and utility for two reasons 
(despite being matter dominated within the concordance model). First, the high redshift 
provides a long lever arm so that slope changes will result in large differences in the distance 
modulus; for example a comparison between the concordance and 'Reiss' cosmologies at 
z — 1.7 has only a difference of 0.15 mag while at z — 6.6 the difference is 1.00 mag. Second, 
it is good science to explore previously unexplored regions and to thoroughly test any recent 
concordance models, especially when the data is arriving even now and for free from Swift. 

(5) I argue that GRBs do not suffer from any evolution effects because the GRB lu- 
minosity indicators are based on light travel times, energy conservation, and the degree of 
relativistic beaming, none of which vary with metallicity or age of the universe. The pop- 
ulation of bursts might evolve in average luminosity, but the luminosity indicators will still 
yield the correct luminosity of each individual burst. 
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(6) I find that the current GRB HD is consistent with the concordance cosmology of 
w = — 1, Q M = 0.27, and £l\ = 0.73. That is, there is no evidence for the Cosmological 
Constant being inconstant. 

This work is only possible from the results of many workers over the past decade and 
more who have spent vast times and efforts accumulating the observations summarized in 
the data tables. These includes the designers, builders, and analysts for the burst detectors 
onboard the GRO, Konus, SAX, HETE-2, Integral, and Swift satellites as well as the many 
optical observers who have sleepless nights to discover the optical transients, measure their 
light curves, and get their redshifts. The key task of this enterprise is the fast distribution 
of burst information via the GCN by Scott Barthelmy. 

I thank Hans Krimm for providing many spectral fit parameters from Swift data, Sergey 
Golenetskii for passing along various E pea k and P values from the Konus- Wind experiment, 
Sheila McBreen for providing the INTEGRAL light curve and fit values for GRB 050502, 
and Hugo Martel for sending numerical results of his gravitational lensing probability dis- 
tributions. I also thank Jochen Greiner for maintaining his web page with an exhaustive 
collection of the data in a convenient format. Neil Gehrels, Hans Krimm, and David Palmer 
helped in improving the manuscript for this paper. 
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Table 1. Prior GRB Hubble Diagrams. 



Reference 


GRBs 


Relations Used a 


Number of Indicators 


Schacfcr (2003a) 


9 


Tlag, V 


18 


Bloom, Frail, & Kulkarni (2003) 


16 b 


E^ 


16 


Xu, Dai, & Liang (2005) 


17 


tjet 


17 


Firmani et al. (2005) 


15 


tjet 


15 


Liang & Zhang (2005) 


15 




15 


Firmani et al. (2006) 


19 


To. 45 


19 


Schaefer (this work) 


69 


Tlag, V, Epeak, tjet, TRT 


243 



a The first five relations in section 3 are indicated by r; a9 , V, E pea k, tjet, and trt respectively. 
The E~/ indicates the tj e t relation but without the later improvement of the correlation with -B pea fe 
The T0.45 (a measure of the burst duration) indicates a relation between the product of powers of 
Epeak an d T3.45 and the isotropic luminosity. 

b An additional 13 GRBs arc discussed but these either do not have redshifts or have only limits 
on the jet break time. 
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Table 2 — Continued 





Experiment 


Z 


tier. 


pb c 


Emin — Emax (keV) 


net. 


»peofe (KeV J 


a 
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050126 
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1.29 


50 


0.70 ± [0.07] 
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51 


47-8 + 23 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


52 
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1.44 


53 
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54 
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99 
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95 
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Swift 
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100 
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101 


237 - 51 + 76 


-1.29 ± 0.07 
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Swift 


4.05 


103 


2.8 ± 0.1 
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Swift 
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105 


2.8 ± 0.2 
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150 
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-1.18 ± 0.23 


[-2.2] ± [0.4] 


95 


060223 


Swift 


4.41 


107 


1.4 ± 0.1 
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150 


108 
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-1.18 ± 0.31 
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95 
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Swift 
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6.7 ± 0.2 
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150 
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Swift 
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Swift 
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95 
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Swift 
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[-2.2] ± [0.4] 


95 
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Swift 
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119 


0.6 ± 0.1 
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150 
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0.5 ± 0.1 
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95 
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125 


1.4 ± 0.1 
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-1.06 ± 0.18 


[-2.2] ± [0.4] 


95 
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Sakamato ct al. 2006b; ( 64 ) Prochaska et al. 2005a; ( 65 ) Gotz ct al. 2005; ( 66 ) McBreen 2005; ( 67 ) Berger et al. 2005b; ( 68 ) Hullinger et al. 2005; 
( 69 ) Foley ct al. 2005; ( 70 ) Cummings ct al. 2005; ( 71 ) Berger & Becker 2005; ( 72 ) Fcnimore et al. 2005a; ( 73 ) Fynbo ct al. 2005c; ( 74 ) Palmer 
ct al. 2005; ( 75 ) Prochaska et al. 2005b; ( 76 ) Cummings et al. 2006a; ( 77 ) Fynbo ct al. 2005d; ( 78 ) Krimm ct al. 2006b; ( 79 ) Crew et al. 2005; 
( 80 ) Kawai ct al. 2005; ( 81 ) Sakamato ct al. 2006a; ( 82 ) Fugazza et al. 2005; ( 83 ) Sato ct al. 2006a; ( 84 ) Jakobsson ct al. 2005a; ( 85 ) Krimm 
ct al. 2006c; ( 86 ) Gal- Yam ct al. 2005; ( 87 ) Golenetskii et al. 2005a; ( 88 ) Quimby et al. 2005a; ( 89 ) Fcnimore ct al. 2005b; ( 90 ) Golenetskii et al. 
2005b; ( 91 ) Hill et al. 2005; ( 92 ) Krimm et al. 2006d; ( 93 ) Melandri et al. 2006; ( 94 ) Sakamato ct al. 2006c; ( 95 ) Krimm 2006; ( 96 ) Piranomontc 
ct al. 2006; ( 97 ) Barbicr et al. 2006; ( 98 ) Grazian ct al. 2006; ( 99 ) Barthelmy 2006a; ( 100 ) Ccnko ct al. 2006; ( 101 ) Fenimore ct al. 2006; ( 102 
) Golenetskii ct al. 2006a; ( 103 ) Fynbo et al. 2006; ( 104 ) Palmer ct al. 2006; ( 105 ) Cucchiara ct al. 2006a; ( 106 ) Sakamoto et al. 2006d; ( 107 ) 
Berger 2006; ( 108 ) Cummings et al. 2006b; ( 109 ) Dupre ct al. 2006; ( 110 ) Cummings ct al. 2006c; ( 111 ) Golenetskii ct al. 2006b; ( 112 ) Cucchiara 
ct al. 2006b; ( 113 ) Parsons et al. 2006a; ( 114 ) Price ct al. 2006; ( 115 ) Barthelmy et al. 2006b; ( 117 ) Berger & Gladders 2006; ( 118 ) Markwardt ct 
al. 2006; ( 119 ) Castro-Tirado ct al. 2006; ( 120 ) Parsons ct al. 2006b; ( 122 ) Peterson et al. 2006; ( 123 ) Sato ct al. 2006b; ( 125 ) Lcdoux ct al. 2006; 
( 126 ) Tueller et al. 2006; ( 127 ) Price et al. 2002a 

b The reported P value is either in units of photons cm -2 s~ 1 for values larger than 0.01 (not in scientific notation) or in units of erg cm -2 s~ 1 for values 
smaller than 0.01 (those in scientific notation) 

c The values reported in square brackets are conservative estimations for uncertainties not reported in the original paper. 

d The values reported in square brackets are average values for a and /3 for cases where no specific value is reported in the literature. 

e The P value for GRB 991208 is reported for a peak time interval of 0.25 seconds, so the peak flux must be multiplied by a factor of 0.85 to convert to 
a flux over the peak 1 second time interval. 

f The P value for GRB 000926 is reported for a peak time interval of 0.25 seconds, so the peak flux must be multiplied by a factor of 0.57 to convert to a 
flux over the peak 1 second time interval. 

s The P value for GRB 040924 is reported for a peak time interval of 0.064 seconds, so the peak flux must be multiplied by a factor of 0.44 to convert to 
a flux over the peak 1 second time interval. 

h The spectral fit parameters were based on a joint fit to the data from Swift and Konus. 

'The P value for GRB 051022 is reported for a peak time interval of 0.064 seconds, so the peak flux must be multiplied by a factor of 0.84 to convert to 
a flux over the peak 1 second time interval. 
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a The uncertainties given in square brackets are conservative estimates for when no error bar is quoted in the original literature. 

b References. — ( 1 ) Jimenez, Band, & Piran 2001; ( 2 ) Bloom ct al. 2003b; ( 3 ) Ghirlanda et al. 2004; ( 4 ) Nava ct al. 2006; ( 5 ) Kulkarni et al. 
1999; ( 6 ) Israel ct al. 1999; ( 7 ) Amati ct al. 2002; ( 8 ) Bjornsson et al. 2001; ( 11 ) Frontera et al. 2001; ( 12 ) Jakobsson ct al. 2003; ( 13 ) Sakamoto 
ct al. 2005a; ( 14 ) Berger et al. 2003; ( 15 ) Price et al. 2003a; ( 16 ) Holland et al. 2003; ( 17 ) Klose et al. 2004; ( 18 ) Andersen et al. 2003; ( 19 
) Jakobsson et al. 2004; ( 20 ) Galassi et al. 2004; ( 21 ) Stanck ct al. 2005; ( 22 ) Krimm ct al. 2005; ( 23 ) Perri et al. 2005; ( 24 ) Hullinger ct al. 
2005; ( 25 ) Hurkctt et al. 2006; ( 26 ) Cummings ct al. 2005; ( 27 ) Blustin ct al. 2006; ( 28 ) Sakamoto et al. 2006a; ( 29 ) Tagliafcrri ct al. 2005; ( 30 
) Frail ct al. 2006; ( 33 ) This work ( 34 ) Golcnetskii et al. 2005; ( 35 ) Racusin ct al. 2005; ( 36 ) Golenetskii et al. 2006; ( 37 ) Romano ct al. 2006; ( 



38 ) Sakamoto et al. 2006d; ( 39 ) Dai & Stanck 2006; ( 40 ) Markwardt ct al. 2006; ( 41 ) Moretti; et al. 2006; 
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Table 4. Luminosities and Luminosity Indicators. 
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13 + 35 


0.5 ± 0.3 


2 


050408 


1.24 


1.1E-6 


± 


2.1E-7 






0.25 ± 0.1 














0.25 ± 0.08 


2 


050416 


0.65 


5.3E-7 


± 


8.5E-8 














15 - 


2.7 - 


f 2.3 


0.51 ± 0.3 


1 


050502 


3.79 


4.3E-7 


± 


1.2E-7 






0.2 ± 0.2 


0.0221 


± 


0.0029 


93 - 


35 + 55 


0.4 ± 0.2 


4 


050505 


4.27 


3.2E-7 


± 


5.4E-8 


6.20E-6 ± 8.5E-7 


0.0014 ± 0.0007 




0.0035 


± 


0.0019 


70 - 


24 ^ 


- 140 


0.40 ± 0.15 


7 


050525 


0.61 


5.2E-6 


± 


7.2E-8 


2.59E-5 ± 1.3E-6 


0.0025 ± 0.0010 


0.11 ± 0.02 


0.0135 


± 


0.0003 


81 - 


1.4 - 


•h 1-4 


0.32 ± 0.03 


2 


050603 


2.82 


9.7E-6 


± 


6.0E-7 






0.03 ± 0.03 


0.0163 


± 


0.0015 


344 


- 52 


+ 52 


0.17 ± 0.02 


1 


050802 


1.71 


5.0E-7 


± 


7.3E-8 








0.0046 


± 


0.0053 








0.8 ± 0.2 


3 


050820 


2.61 


3.3E-7 


± 


5.2E-8 






0.7 ± 0.3 








246 


- 40 


+ 76 


2.0 ± 0.5 




050824 


0.83 


9.3E-8 


± 


3.8E-8 




















11 ± 2 


1 


050904 


6.29 


2.5E-7 


± 


3.5E-8 


2.00E-5 ± 2.0E-6 


0.0097 ± 0.0024 




0.0023 


± 


0.0026 


436 


- 90 


+ 200 


0.6 ± 0.2 


7 


050908 


3.35 


9.8E-8 


± 


1.5E-8 














41 - 


5 + 


9 


1.5 ± 0.3 


1 


050922 


2.20 


2.0E-6 


± 


7.3E-8 






0.06 ± 0.02 


0.0033 


± 


0.0006 


198 


- 22 


+ 38 


0.13 ± 0.02 


2 


051022 


0.80 


1.1E-5 


± 


8.7E-7 


3.40E-4 ± 1.2E-5 


0.0029 ± 0.0001 




0.0122 


± 


0.0004 


510 


- 20 


+ 22 


0.19 ± 0.04 


9 


051109 


2.35 


7.8E-7 


± 


9.7E-8 














161 


- 35 


+ 130 


1.3 ± 0.4 


3 


051111 


1.55 


3.9E-7 


± 


5.8E-8 






1.02 ± 0.1 


0.0024 


± 


0.0007 








3.2 ± 1 


1 


060108 


2.03 


1.1E-7 


± 


1.1E-7 








0.0032 


± 


0.0058 


65 - 


10 ^ 


- 600 


0.4 ± 0.2 


2 


060115 


3.53 


1.3E-7 


± 


1.6E-8 














62 - 


6 + 


19 


0.4 ± 0.2 


4 


060116 


6.60 


2.0E-7 


± 


1.1E-7 














139 


- 36 


+ 400 


1.3 ± 0.5 


9 


060124 


2.30 


1.1E-6 


± 


1.2E-7 


3.37E-5 ± 3.4E-6 


0.0021 ± 0.0002 


0.08 ± 0.04 


0.0140 


± 


0.0020 


237 


- 51 


+ 76 


0.3 ± 0.1 


4 


060206 


4.05 


4.4E-7 


± 


1.9E-8 






0.1 ± 0.1 


0.0025 


± 


0.0016 


75 - 


12 ^ 


- 12 


1.25 ± 0.25 


1 


060210 


3.91 


5.5E-7 


± 


2.2E-8 


1.94E-5 ± 1.2E-6 


0.0005 ± 0.0001 


0.13 ± 0.08 


0.0019 


± 


0.0004 


149 


- 35 + 400 


0.5 ± 0.2 


9 


060223 


4.41 


2.1E-7 


± 


3.7E-8 






0.38 ± 0.1 


0.0075 


± 


0.0033 


71 - 


io ^ 


- 100 


0.50 ± 0.1 


4 


060418 


1.49 


1.5E-6 


± 


5.9E-8 






0.26 ± 0.06 


0.0070 


± 


0.0005 


230 


- [20] + [20] 


0.32 ± 0.08 


3 


060502 


1.51 


3.7E-7 


± 


1.6E-7 






3.5 ± 0.5 


0.0010 


± 


0.0017 


156 


- 33 


+ 400 


3.1 ± 0.3 


1 


060510 


4.90 


1.0E-7 


± 


1.7E-8 








0.0028 


± 


0.0019 


95 - 


[30] 


+ [60] 




9 


060526 


3.21 


2.4E-7 


± 


3.3E-8 


1.17E-6 ± 1.7E-7 


0.0034 ± 0.0014 


0.13 ± 0.03 


0.0112 


± 


0.0039 


25 - 


[5] - 


- [5] 


0.20 ± 0.05 


3 


060604 


2.68 


9.0E-8 


± 


1.6E-8 






5 ± 1 








40 - 


[5] - 


- [5] 


0.6 ± 0.2 


2 


060605 


3.80 


1.2E-7 


± 


5.5E-8 






5 ± 3 








169 


- [30] + [200] 


2.0 ± 0.5 


2 


060607 


3.08 


2.7E-7 


± 


8.1E-8 






2.0 ± 0.5 


0.0059 


± 


0.0014 


120 


- 17 


+ 190 


2.0 ± 0.2 


1 



a The uncertainties given in square brackets are conservative estimates for when no error bar is quoted in the original literature. 
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Table 5. Light Curve Properties. 



Experiment For V, trt, ^peak F° r T lag Time Resolution (sec) 

BATSE 25-320 keV 25-50 and 100-320 keV 0.064 

Konus 50-200 keV only one channel available 0.064 

HETE-2 30-400 keV 7-40 and 30-400 keV 0.164 

Swift 25-350 keV 25-50 and 100-350 keV 0.064 
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Table 6. Derived Distance Moduli. 



GRB 


z 


Ml a 




M2 a 








M3 a 








M4 a 








M5 a 








a 






H 


b 




970228 


0.70 




42. 


.OO 


1 

± 


1 


31 


42 


. 1 1 


± 


1 


.16 












43 


1 O 

.42 


± 


1 


.21 


42.79 


zb 


0.71 


42 


56 


zb 


0.70 


970508 


0.84 


43.09 ± 1.19 


43. 


"1 Q 

. lo 


1 

zt 


1 


.34 


45 


.00 


zt 


1 


.04 


43 


.\)Z 


1 

zt 





.63 


43 


no 

,uy 


1 

zt 


1 


.20 


43.92 


zb 


0.43 


43 


68 


zb 


0.44 


970828 


0.96 




42. 


on 


1 

Zt 


1 


28 


44 


C\ r 7 


1 

zt 


1 


.03 


43 


.01 


1 

± 





.63 


43 


07 
,z 1 


zt 


1 


.25 


43.58 


zb 


0.46 


43 


.32 


zb 


0.48 


971214 


3.42 


48.76 zb 1.50 


48. 


fin 
Oy 


1 

Zt 


1 


.19 


47 


r, it 


zt 


1 


.04 












49 


OO 


zt 


1 


.33 


48.44 


zb 


0.62 


48 


06 


zb 


0.61 


980613 


1.10 














45 


Q It 
.OO 


zt 


1 


.33 






















45.85 


zb 


1.33 


45 


.63 


zb 


1.30 


980703 


0.97 


44.47 ± 1.03 


44. 


(VI 

yi 


1 

Zt 


1 


.28 


46 


t~\l-. 

.UO 


zt 


1 


.02 


43 


",n 

.Oy 


1 

± 





.58 


42 


1 "1 
.11 


zt 


1 


.20 


44.07 


zb 


0.40 


43 


91 


zb 


0.42 


990123 


1.61 


43.23 ± 1.02 


44. 


OO 


1 

Zt 


1 


.21 


45 


~,n 

,0y 


zt 


1 


.05 


45 


.04 


1 

± 





66 












44.99 


zb 


0.45 


44. 


.63 


zb 


0.47 


990506 


1.31 


44.82 ± 1.15 


44. 


O -1 

.24 


1 

Zt 


1 


.23 


44 


1 n 

.ly 


zt 


1 


.03 












43. 


,yO 


zt 


1 


.22 


44.31 


zb 


0.57 


43. 


.99 


zb 


0.57 


990510 


1.62 


46.59 ± 1.08 


45. 


07 
Z i 


_l_ 
Zt 


1 


23 


44 


0f[ 
.ZD 


_l_ 

zt 


1 


.01 


45 




_l_ 

zt 





.56 


45 


KO 


_l_ 

zt 


1 


.22 


45.46 


zb 


0.40 


45. 


.20 


zb 


0.41 


990705 


0.84 




45. 


K 
ZD 


_l_ 
zt 


1 


.22 


43 


.■Jo 


zt 


1 


.01 


12 


on 


_|_ 

zt 





62 


45 


.ol 


zt 


1 


.32 


43.77 


zb 


0.46 


43. 


.59 


zb 


0.47 


990712 


0.43 














11 


QQ 
■ OO 


zt 


1 


.03 


11 


.04 


1 

zt 





62 












41.63 


zb 


0.53 


41. 


.50 


zb 


0.56 


991208 


0.71 




40. 


,0U 


zt 


1 


.32 


42 


1 Q 
.lo 


zt 


1 


.02 












42 


n 1 

Ul 


zt 


1 


.21 


41.69 


zb 


0.67 


41 


46 


zb 


0.67 


991216 


1.02 


43.55 ± 1.08 


42. 


CO 

OO 


1 

zt 


1 


.23 


42 


■ I Z 


zt 


1 


.03 


43 


.OO 


1 

zt 





.64 


43. 


,04 


zt 


1 


.25 


43.32 


zb 


0.43 


43. 


.00 


zb 


0.44 


000131 


4.50 




47. 


i n 
1U 


1 

zt 


1 


.25 


47 


■ (\J 


zt 


1 


.04 












48 


An 

,4y 


zt 


1 


.38 


47.72 


zb 


0.69 


47. 


.37 


zb 


0.68 


000210 


0.85 




40. 


on 


1 

zt 


1 


.32 


43 




1 

zt 


1 


.02 












41 


no 

,yo 


zt 


1 


.21 


42.45 


zb 


0.67 


42. 


18 


zb 


0.67 


000911 


1.06 




44 


00 


1 

zt 


1 


21 


45 


.0 i 


1 

zt 


1 


.06 












44 


QA 

,o4 


1 

zt 


1 


.32 


45.10 


zb 


0.68 


44. 


70 


zb 


0.68 


000926 


2.07 




46. 


Zo 


_l_ 
zt 


1 


23 


14 


on 


zt 


1 


.02 












47. 


Q O 

oy 


zt 


1 


.45 


45.62 


zb 


0.69 


45 


.32 


zb 


0.68 


010222 


1.48 




43. 


OO 


zt 


1 


.23 


43. 


.O ( 


1 

zt 


1 


.02 


45 


.uo 


_l_ 
zt 





.53 


43 


. { 1 


zt 


1 


.25 


44.48 


zb 


0.42 


44 


10 


zb 


0.45 


010921 


0.45 


42.81 ± 1.06 


40. 


(17 

y r 


1 

zt 


1 


.98 


43. 


"1 It 
.lO 


zt 


1 


.03 












41 


"i 7 
. 1 ( 


zt 


1 


.20 


42.35 


zb 


0.60 


42. 


.22 


zb 


0.59 


011211 


2.14 














17 


ri7 

,U f 


zt 


1 


.03 


45 


. 1 1 


zt 





63 












45.64 


zb 


0.53 


45 


.52 


zb 


0.56 


020124 


3.20 


47.83 zt 1.23 


48. 


51 


1 


1 


.26 


16 


.25 


1 

± 


1 


.05 


16 


.53 


1 





64 


47. 


.36 


1 

± 


1 


.24 


46.98 


zb 


0.44 


46. 


.75 


zb 


0.45 


020405 


0.70 




44. 


AO 

Uo 


1 

zt 


1 


25 


44 


,0U 


1 

zt 


1 


.11 


43 


.04 


1 

zt 





.77 


42 


an 

oy 


1 

zt 


1 


.21 


43.68 


zb 


0.51 


43 


43 


zb 


0.52 


020813 


1.25 


44.39 ± 1.04 


45. 


oz 


1 

zt 


1 


23 


14 


no 


zt 


1 


.02 


14 


O 1 


1 

zt 





,60 


42. 


no 


zt 


1 


.20 


44.11 


zb 


0.41 


43. 


.89 


zb 


0.42 


020903 


0.25 














10 


7fl 
. {)) 


zt 


1 


.14 






















40.79 


zb 


1.14 


40 


.93 


zb 


1.13 


021004 


2.32 


46.43 zt 1.26 


46. 


7^ 
1 


zt 


2 


.06 


46 


. ( 


zt 


1 


.16 


45 


,o4 


zt 





.81 


47. 


■ On 


zt 


1 


.34 


46.45 


zb 


0.52 


46. 


.26 


zb 


0.52 


021211 


1.01 


44.08 ± 1.00 












42 


.32 


1 

± 


1 


.02 












44 


.62 


1 

± 


1 


.21 


43.57 


zb 


0.61 


43 


.37 


zb 


0.61 


030115 


2.50 


46.58 ± 1.15 


47. 


-1 I 


1 

zt 


1 


.61 


46 


. 04 


zt 


1 


.13 












45 


.01 


1 

zt 


1 


.28 


46.44 


zb 


0.63 


46. 


21 


zb 


0.62 


030226 


1.98 


46.93 zt 1.50 


47. 


z 1 


zt 


1 


.71 


46 


. ( O 


zt 


1 


.07 


46 


OA 


zt 





.66 


46 


AQ 

/ly 


_l_ 

zt 


1 


.26 


46.51 


zb 


0.47 


46. 


.32 


zb 


0.48 


030323 


3.37 














46 


.86 


zb 


1 


.56 












47. 


.37 


zb 


1 


.46 


47.13 


zb 


1.07 


46. 


.90 


zb 


1.05 


030328 


1.52 


45.25 ± 1.49 


44. 


77 


± 


1 


.29 


44 


.99 


± 


1 


.02 


44 


.63 


zb 





.60 












44.78 


zb 


0.46 


44. 


.58 


zb 


0.48 


030329 


0.17 


42.04 zt 1.04 


41. 


70 


± 


1 


28 


39 


.70 


zb 





.98 


38 


.92 


zb 





.53 


40 


,64 


zb 


1 


.20 


39.88 


zb 


0.38 


39. 


.83 


zb 


0.40 


030429 


2.66 




47. 


.81 


± 


1 


.88 


45 


.57 


zb 


1 


.11 


46 


.57 


zb 





.84 


46 


.73 


zb 


1 


.26 


46.46 


zb 


0.56 


46. 


.27 


zb 


0.57 


030528 


0.78 


42.80 zt 1.02 


44 


90 


± 


1 


.82 


14 


.33 


zb 


1 


.04 












46 


.23 


zb 


1 


.25 


44.30 


zb 


0.59 


44. 


.22 


zb 


0.59 


040924 


0.86 


44.02 zt 1.01 












42 


.72 


zb 


1 


.01 












45 


.22 


zb 


1 


.21 


43.84 


zb 


0.61 


43. 


.62 


zb 


0.61 


041006 


0.71 


45.08 zt 1.26 


44. 


.24 


± 


1 


.27 


42 


.51 


zb 


1 


.06 


14 


.21 


zb 





.68 


43 


.41 


zb 


1 


.23 


43.91 


zb 


0.45 


43 


72 


zb 


0.46 


050126 


1.29 


45.48 zt 1.01 


46. 


92 


± 


1 


.44 


45 


.93 


zb 


1 


.04 












44 


,87 


zb 


1 


.22 


45.71 


zb 


0.57 


45. 


.58 


zb 


0.57 


050318 


1.44 




46. 


46 


± 


1 


.28 


44 


.31 


zb 


1 


.04 


45 


.91 


zb 





.67 


46. 


.27 


zb 


1 


.21 


45.71 


zb 


0.47 


45. 


.51 


zb 


0.49 


050319 


3.24 




46. 


64 


± 


1 


.71 






















48. 


.82 


zb 


1 


.25 


48.06 


zb 


1.01 


47. 


73 


zb 


0.99 


050401 


2.90 


46.20 ± 1.21 


47. 


12 


± 


1 


.22 


45 


.36 


zb 


1 


.05 












48. 


,74 


zb 


1 


.31 


46.67 


zb 


0.59 


46. 


.32 


zb 


0.58 


050406 


2.44 


48.22 ± 1.24 












46 


.48 


zb 


1 


.39 












49 


.09 


zb 


1 


.45 


47.92 


zb 


0.78 


47. 


.71 


zb 


0.76 


050408 


1.24 


45.23 ± 1.10 
































45. 


.88 


zb 


1 


.28 


45.51 


zb 


0.84 


45 


.25 


zb 


0.83 


050416 


0.65 














41 


.51 


zb 


1 


.04 












45 


36 


zb 


1 


.43 


42.84 


zb 


0.84 


42 


76 


zb 


0.83 


050502 


3.79 


47.35 zt 1.52 


50. 


15 


± 


1 


.23 


46 


.98 


zb 


1 


.26 












47. 


30 


zb 


1 


.40 


48.06 


zb 


0.67 


47. 


.72 


zb 


0.66 


050505 


4.27 




47. 


11 


± 


1 


.52 


46 


.97 


zb 


1 


.20 


48 


.57 


zb 





.97 


47. 


,74 


zb 


1 


.31 


47.76 


zb 


0.60 


47 


,48 


zb 


0.60 


050525 


0.61 


44.10 ± 1.02 


44. 


40 


± 


1 


.24 


42 


.05 


zb 





.99 


43 


.26 


zb 





.66 


43 


.45 


zb 


1 


.20 


43.34 


zb 


0.42 


43 


,14 


zb 


0.43 


050603 


2.82 


45.81 ± 1.50 


45. 


75 


± 


1 


.21 


45 


.58 


zb 


1 


.06 












44. 


,74 


zb 


1 


.21 


45.45 


zb 


0.61 


45. 


.08 


zb 


0.60 
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Table 6 — Continued 



GRB 


2 












M2 a 






M3 a 






M4 a 


M5 a 






H 


a 




H 


b 




050802 


1.71 












45. 


87 


± 


1 


,96 












45.47 


± 


1 


.24 


45.59 


± 


1.05 


45.38 


± 


1.03 


050820 


2.61 


45. 


.97 


± 


1 


.11 












48.56 


± 


1 


.07 




45.12 


± 


1 


.24 


46.70 


± 


0.66 


46.45 


± 


0.65 


050824 


0.83 
































43.35 


± 


1 


.29 


43.35 


± 


1.29 


43.35 


± 


1.29 


050904 


6.29 












47. 


18 


± 


1 


,93 


51.15 


± 


1 


.12 


49.36 ± 0.72 


47.89 


± 


1 


.28 


49.34 


± 


0.53 


49.00 


± 


0.55 


050908 


3.35 






















46.92 


± 


1 


.03 




47.05 


± 


1 


.23 


46.97 


± 


0.79 


46.76 


± 


0.79 


050922 


2.20 


46. 


.57 


± 


1 


.07 


44. 


05 


± 


1 


.32 


45.97 


± 


1 


.03 




46.58 


± 


1 


.22 


45.92 


± 


0.57 


45.60 


± 


0.57 


051022 


0.80 












43, 


62 


± 


1 


.25 


44.80 


± 


1 


.03 


43.88 ± 0.53 


43.48 


± 


1 


.23 


43.95 


± 


0.42 


43.69 


± 


0.45 


051109 


2.35 






















46.70 


± 


1 


.09 




44.64 


± 


1 


.26 


45.82 


± 


0.83 


45.56 


± 


0.82 


051111 


1.55 


14. 


96 


± 


1 


.01 


44. 


.77 


± 


1 


,42 












43.84 


± 


1 


.26 


44.58 


± 


0.69 


44.44 


± 


0.68 


060108 


2.03 












47. 


04 


± 


2 


,61 


47.00 


± 


1 


.50 




48.18 


± 


1 


.74 


47.43 


± 


1.04 


47.19 


± 


1.03 


060115 


3.53 






















47.42 


± 


1 


.02 




48.49 


± 


1 


.37 


47.80 


± 


0.82 


47.54 


± 


0.81 


060116 


6.60 






















49.40 


± 


1 


.27 




47.19 


± 


1 


.42 


48.42 


± 


0.95 


48.10 


± 


0.94 


060124 


2.30 


46. 


96 


± 


1 


.15 


47. 


.57 


± 


1 


,24 


47.03 


± 


1 


.10 


47.11 ± 0.66 


46.20 


± 


1 


.28 


47.02 


± 


0.44 


46.73 


± 


0.45 


060206 


4.05 


48. 


16 


± 


1 


.49 


46, 


04 


± 


1 


,59 


46.69 


± 


1 


.05 




45.86 


± 


1 


.22 


46.63 


± 


0.64 


46.37 


± 


0.63 


060210 


3.91 


47. 


59 


± 


1 


.21 


45. 


21 


± 


1 


.34 


47.63 


± 


1 


.10 


49.58 ± 0.70 


46.78 


± 


1 


.31 


48.08 


± 


0.46 


47.70 


± 


0.47 


060223 


4.41 


47. 


56 


± 


1 


.05 


49, 


09 


± 


1 


.41 


47.51 


± 


1 


.28 




47.95 


± 


1 


.24 


47.93 


± 


0.61 


47.65 


± 


0.60 


060418 


1.49 


44. 


97 


± 


1 


.03 


45, 


31 


± 


1 


.26 


46.08 


± 


1 


.02 




45.36 


± 


1 


.24 


45.45 


± 


0.56 


45.18 


± 


0.56 


060502 


1.51 


43. 


66 


± 


1 


.10 


43, 


12 


± 


2 


,41 


46.93 


± 


1 


.17 




43.93 


± 


1 


.28 


44.71 


± 


0.66 


44.55 


± 


0.65 


060510 


4.90 












48, 


,17 


± 


1 


,62 


49.00 


± 


1 


.18 












48.71 


± 


0.95 


48.42 


± 


0.93 


060526 


3.21 


48. 


32 


± 


1 


.05 


49, 


23 


± 


1 


,34 


44.99 


± 


1 


.06 


46.94 ± 0.77 


48.68 


± 


1 


.25 


47.35 


± 


0.46 


47.13 


± 


0.47 


060604 


2.68 


45. 


22 


± 


1 


.03 












46.69 


± 


1 


.03 




48.12 


± 


1 


.28 


46.48 


± 


0.63 


46.31 


± 


0.63 


060605 


3.80 


45. 


16 


± 


1 


.28 












49.44 


± 


1 


.17 




46.54 


± 


1 


.32 


47.21 


± 


0.72 


46.97 


± 


0.72 


060607 


3.08 


45. 


17 


± 


1 


.08 


47, 


83 


± 


1 


,34 


47.70 


± 


1 


.09 




45.51 


± 


1 


.24 


46.48 


± 


0.59 


46.26 


± 


0.58 



a This assumes the concordance cosmology with Qm = 0.27 and ui = —1. The units of the distance moduli are magnitudes. 
b This assumes the best fit cosmology from the gold sample of Riess et al. (2004) with wq = —1.31 and w' = 1.48. 
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Table 7. Gravitational Lensing and Malmquist Biases. 



GRB 


Experiment 




z 




L (erg 


s 


-l 


Pbolo 


(erg cm 2 s - 1 ) 


fJ-bias (mag) 


970228 


Konus 





.695 


1 


7E+52 


± 


4 


.OE+51 


7. 


3E- 


■6 


± 


4. 


.3E-7 


0.05 ± 0.01 


970508 


BATSE 





.835 


6 


5E+51 


± 


9 


.5E+50 


3 


3E- 


■6 


± 


3 


.3E-7 


0.01 ± 0.01 


970828 


Konus 





.9578 


8. 


5E+51 


± 


1 


.8E+51 


1 


0E- 


■5 


± 


1 


1E-6 


0.05 ± 0.01 


971214 


BATSE 


3 


.418 


2 


1E+53 


± 


6 


.4E+52 


7. 


.5E- 


■7 


± 


2 


.4E-8 


-0.15 ± 0.02 


980613 


SAX 


1 


.0964 


8. 


1E+51 


± 


7 


.OE+51 


3 


0E- 


■7 


± 


8. 


.3E-8 


0.08 ± 0.03 


980703 


BATSE 





.966 


2 


.7E+51 


± 


3 


.9E+50 


1 


.2E- 


6 


± 


3 


6E-8 


-0.01 ± 0.01 


990123 


BATSE 


1 


.61 


3 


OE+52 


± 


5 


.7E+51 


1 


3E- 


■5 


± 


5 


.OE-7 


0.06 ± 0.01 


990506 


BATSE 


1 


.30658 


6 


1E+52 


± 


1 


.3E+52 


1 


1E- 


■5 


± 


1 


.5E-7 


0.07 ± 0.01 


990510 


BATSE 


1 


.619 


1 


6E+52 


± 


2 


.6E+51 


3 


3E- 


■6 


± 


1 


.2E-7 


0.01 ± 0.01 


990705 


Konus 





.8424 


5. 


8E+51 


± 


1 


.3E+51 


6 


6E- 


6 


± 


2 


.6E-7 


0.03 ± 0.00 


990712 


SAX 





.43 


2. 


3E+50 


± 


6 


.6E+49 


3 


.5E- 


■6 


± 


2 


.9E-7 


0.04 ± 0.01 


991208 


Konus 





.707 


1 


5E+52 


± 


3 


.OE+51 


2 


1E- 


■5 


± 


2 


1E-6 


0.05 ± 0.01 


991216 


BATSE 


1 


.02 


2 


.4E+52 


± 


4 


.6E+51 


4. 


1E- 


■5 


± 


3 


.8E-7 


0.04 ± 0.01 


000131 


Konus 


4 


.5 


1 


1E+53 


± 


3 


.8E+52 


7. 


3E- 


■7 


± 


8 


.3E-8 


-0.22 ± 0.04 


000210 


Konus 





.846 


2. 


OE+52 


± 


4 


.2E+51 


2 


0E- 


■5 


± 


2 


1E-6 


0.06 ± 0.01 


000911 


Konus 


1 


.0585 


2. 


4E+53 


± 


7 


.9E+52 


1 


9E- 


■5 


± 


1 


9E-6 


0.16 ± 0.03 


000926 


Konus 


2 


.066 


4. 


6E+52 


± 


1 


.5E+52 


2 


9E- 


■6 


± 


2 


.9E-7 


0.03 ± 0.01 


010222 


Konus 


1 


.476 


1 


6E+52 


± 


3 


.2E+51 


2 


3E- 


■5 


± 


7 


.2E-7 


0.04 ± 0.01 


010921 


HETE 





.45 


1 


2E+51 


± 


2 


.OE+50 


1 


.8E- 


6 


± 


1 


.6E-7 


0.03 ± 0.01 


011211 


SAX 


2 


.14 


7. 


0E+50 


± 


1 


.9E+50 


9. 


.2E- 


■8 


± 


9 


3E-9 


-0.07 ± 0.01 


020124 


HETE 


3 


198 


1 


5E+52 


± 


3 


.1E+51 


6 


1E- 


■7 


± 


1 


.OE-7 


0.09 ± 0.01 


020405 


Konus 





.695 


1 


OE+52 


± 


2 


.1E+51 


7. 


.4E- 


6 


± 


3 


.1E-7 


0.04 ± 0.01 


020813 


HETE 


1 


.254 


6 


6E+51 


± 


9 


.4E+50 


3 


.8E- 


■6 


± 


2 


.6E-7 


0.04 ± 0.01 


020903 


HETE 





,25 


8. 


2E+48 


± 


4 


.9E+48 


3 


.4E- 


■8 


± 


8 


8E-9 


-0.76 ± 0.27 


021004 


HETE 


2 


.323 


5. 


.5E+51 


± 


1 


.6E+51 


2 


3E- 


■7 


± 


5 


.5E-8 


0.00 ± 0.01 


021211 


HETE 


1 


.006 


1 


1E+52 


± 


1 


.7E+51 


2 


3E- 


■6 


± 


1 


.7E-7 


0.05 ± 0.01 


030115 


HETE 


2 


.5 


1 


1E+52 


± 


3 


.2E+51 


3 


2E- 


■7 


± 


5 


.1E-8 


0.03 ± 0.01 


030226 


HETE 


1. 


.984 


4. 


2E+51 


± 


9 


.7E+50 


2 


6E- 


■7 


± 


4 


.7E-8 


0.02 ± 0.01 


030323 


HETE 


3 


,372 


1 


1E+52 


± 


6 


.1E+51 


1 


.2E- 


■7 


± 


6 


.OE-8 


-0.33 ± 0.11 


030328 


HETE 


1 


.52 


2. 


3E+51 


± 


5 


.7E+50 


1 


6E- 


■6 


± 


1 


IE- 7 


0.08 ± 0.01 


030329 


HETE 





.1685 


1 


5E+51 


± 


1 


.9E+50 


2 


0E- 


■5 


± 


1 


.OE-6 


0.01 ± 0.01 


030429 


HETE 


2 


.6564 


2 


8E+51 


± 


6 


.5E+50 


2 


0E- 


■7 


± 


5 


.4E-8 


0.00 ± 0.01 


030528 


HETE 





.782 


5. 


0E+50 


± 


7 


.2E+49 


1 


6E- 


■7 


± 


3 


2E-8 


-0.03 ± 0.01 


040924 


HETE 





.859 


1 


7E+52 


± 


2 


.5E+51 


2. 


6E- 


6 


± 


2 


.8E-7 


0.07 ± 0.01 


041006 


HETE 





.712 


2 


3E+51 


± 


4 


.9E+50 


2. 


.5E- 


6 


± 


1 


.4E-7 


0.05 ± 0.01 


050126 


Swift 


1 


.29 


1 


8E+51 


± 


2 


.7E+50 


1 


.1E- 


■7 


± 


1 


.3E-8 


0.03 ± 0.01 


050318 


Swift 


1 


.44 


4. 


5E+51 


± 


8 


.2E+50 


5. 


2E- 


■7 


± 


6 


.3E-8 


0.08 ± 0.01 


050319 


Swift 


3 


.24 


7. 


9E+52 


± 


2 


.6E+52 


2 


3E- 


■7 


± 


3 


.6E-8 


-0.02 ± 0.01 


050401 


Swift 


2 


.9 


1 


4E+53 


± 


3 


.9E+52 


2. 


1E- 


■6 


± 


2 


.2E-7 


0.13 ± 0.03 


050406 


Swift 


2 


.44 


9. 


3E+51 


± 


4 


.OE+51 


4. 


2E- 


■8 


± 


1 


1E-8 


-0.20 ± 0.05 


050408 


HETE 


1 


.236 


2. 


2E+52 


± 


6 


.9E+51 


1 


.IE- 


6 


± 


2 


.1E-7 


0.16 ± 0.03 


050416 


Swift 





.6535 


6 


9E+50 


± 


2 


.7E+50 


5. 


3E- 


■7 


± 


8 


.5E-8 


0.15 ± 0.03 


050502 


Integral 


3 


.793 


1 


0E+53 


± 


3 


.7E+52 


4. 


3E- 


■7 


± 


1 


2E-7 


-0.02 ± 0.01 


050505 


Swift 


4 


.27 


1 


2E+52 


± 


3 


.8E+51 


3 


2E- 


■7 


± 


5 


.4E-8 


0.17 ± 0.04 


050525 


Swift 





.606 


4. 


9E+51 


± 


6 


.2E+50 


5. 


.2E- 


6 


± 


7 


.2E-8 


0.03 ± 0.01 


050603 


Swift 


2 


.821 


1 


2E+53 


± 


2 


.8E+52 


9. 


7E- 


■6 


± 


6 


.OE-7 


0.16 ± 0.03 
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Table 7— Continued 



GRB 


Experiment 




z 




L (erg 


s 


-l 


Pbolo (erg cm 


2 s- 1 ) 


ix bias (mag) 


050802 


Swift 


1 


71 


9 


5E+51 


± 


3 


1E+51 


5.0E-7 


± 


7. 


3E- 


•8 


0.19 ± 0.04 


050820 


Swift 


2. 


.612 


1 


OE+52 


± 


2 


5E+51 


3.3E-7 


± 


5. 


2E- 


■8 


0.14 ± 0.02 


050824 


Swift 





.83 


2 


4E+50 


± 


6 


9E+49 


9.3E-8 


± 


3 


8E- 


•8 


0.00 ± 0.01 


050904 


Swift 


6 


29 


6 


OE+52 


± 


1 


7E+52 


2.5E-7 


± 


3 


.5E- 


■8 


0.03 ± 0.01 


050908 


Swift 


3 


.35 


7. 


5E+51 


± 


1 


.8E+51 


9.8E-8 


± 


1 


.5E- 


■8 


0.00 ± 0.01 


050922 


Swift 


2 


.198 


7. 


8E+52 


± 


1 


.6E+52 


2.0E-6 


± 


7. 


.3E- 


■8 


0.09 ± 0.01 


051022 


Konus 


0. 


8 


1 


5E+52 


± 


3 


OE+51 


1.1E-5 


± 


8. 


.7E- 


■7 


0.05 ± 0.01 


051109 


Swift 


2. 


.346 


1 


3E+52 


± 


4. 


1E+51 


7.8E-7 


± 


9 


.7E- 


■8 


0.19 ± 0.04 


051111 


Swift 


1 


.55 


3 


9E+51 


± 


6 


1E+50 


3.9E-7 


± 


5 


.8E- 


•8 


0.09 ± 0.01 


060108 


Swift 


2. 


03 


1 


2E+52 


± 


4. 


.6E+51 


1.1E-7 


± 


1 


.1E- 


■7 


-0.05 ± 0.02 


060115 


Swift 


3 


53 


2. 


1E+52 


± 


7. 


.6E+51 


1.3E-7 


± 


1 


.6E- 


■8 


0.01 ± 0.01 


060116 


Swift 


6 


.6 


3 


9E+52 


± 


1 


.5E+52 


2.0E-7 


± 


1 


.1E- 


■7 


0.03 ± 0.01 


060124 


Swift 


2. 


.297 


2. 


7E+52 


± 


6 


.3E+51 


1.1E-6 


± 


1 


2E- 


■7 


0.12 ± 0.02 


060206 


Swift 


4. 


045 


1 


6E+52 


± 


3 


.7E+51 


4.4E-7 


± 


1 


.9E- 


■8 


0.13 ± 0.02 


060210 


Swift 


3 


91 


1 


6E+52 


± 


4. 


.2E+51 


5.5E-7 


± 


2 


2E- 


•8 


0.16 ± 0.03 


060223 


Swift 


4. 


.41 


3 


4E+52 


± 


7. 


1E+51 


2+E-7 


± 


3 


.7E- 


■8 


0.04 ± 0.01 


060418 


Swift 


1 


49 


2. 


4E+52 


± 


4 


.7E+51 


1.5E-6 


± 


5 


.9E- 


■8 


0.07 ± 0.01 


060502 


Swift 


1 


.51 


1 


7E+51 


± 


2 


2E+50 


3.7E-7 


± 


1 


.6E- 


■7 


0.22 ± 0.05 


060510 


Swift 


4. 


.9 


4. 


OE+52 


± 


2 


4E+52 


1.0E-7 


± 


1 


.7E- 


■8 


-0.29 ± 0.06 


060526 


Swift 


3 


21 


1 


1E+52 


± 


2 


.1E+51 


2.4E-7 


± 


3 


.3E- 


■8 


0.08 ± 0.01 


060604 


Swift 


2. 


.68 


2 


8E+51 


± 


5. 


.5E+50 


9.0E-8 


± 


1 


.6E- 


■8 


0.03 ± 0.01 


060605 


Swift 


3 


.8 


9 


1E+51 


± 


2 


.4E+51 


1.2E-7 


± 


5 


.5E- 


■8 


0.03 ± 0.01 


060607 


Swift 


3 


.082 


6 


OE+51 


± 


8. 


9E+50 


2.7E-7 


± 


8. 


.1E- 


■8 


0.16 ± 0.03 



Table 8. How Robust are the x 2 Contours? 



Case GRBs Best \ 2 Position of long axis Constraint on wo 

x 2 (w = -1) x 2 ('Riess') 



Standard 


69 


70.7 


w a 


+ 2.1TO0 




-1.0 ±0.3 


wo 


= 0.2 ± 0.6 


72.3 


80.3 


Minus 2 highest-2 GRBs 


67 


68.4 


W a 


+ 2.1iu 




-0.8 ±0.3 


wo 


= 0.5 ± 2 


69.5 


73.0 


Plus 5 outliers 


72 


209.9 


W a 


+ 2.2w 




-0.6 ±0.3 


wo 


> 2 


240.9 


247.0 


Plus 4 outliers (without GRB 980425) 


71 


96.7 


W a 


+ 2.2w 




-0.4 ±0.3 


wo 


> 


105.1 


110.2 


Minus GRBs 050802, 060108, 060108 


66 


67.4 


W a 


+ 2.1wj 




-0.8 ±0.3 


WO 


= 0.6 ±0.9 


69.6 


77.7 


Minus Ti a g — L relation 


69 


67.8 


W a 


+ 2.1u;o 




-1.4 ±0.4 


w 


= -0.3 ± 1.2 


67.9 


76.7 


Minus V — L relation 


69 


58.3 


W a 


+ 2.lw 




-0.6 ±0.3 


w 


= 0.2 ± 1.4 


60.3 


64.6 


Minus E pea k — L relation 


67 


86.0 


W a 


+ 2.2w 




-0.4 ±0.3 


wo 


= 0.4 ± 1.3 


88.3 


90.7 


Minus E pea k ~ E-f relation 


69 


75.5 


W a 


+ 2.1w 




-0.4 ±0.3 


w 


= 1.2 ± 1.4 


79.2 


82.4 


Minus tut — L relation 


68 


59.2 


W a 


+ 2.1w 




-0.8 ±0.4 


wo 


= 0.8 ± 1.1 


61.5 


68.8 


Minus V — L and trt — L relations 


66 


51.4 


W a 


+ 2.0w 




-0.6 ±0.3 


Wo 


= 0.8 ± 1.2 


54.0 


57.7 


Ti a g — L relation only 


39 


25.2 


W a 


+ 2.Qw 




2.5 ±0.7 


w a 


> 2 


37.1 


33.0 


V — L relation only 


51 


48.0 


W a 


+ 2.lw 




-0.6 ±0.7 


wo 


> 


49.7 


50.2 


Epeak ~ L relation only 


62 


60.8 


W a 


+ 2.lw 




0.0 ±0.4 


Wo 


= 1.6 ±2.0 


63.4 


63.3 


Epeak — E~{ relation only 


27 


25.0 


W a 


+ 2.0w 




-3 ±3 


w 


= -0.5 ±2.5 


25.0 


27.6 


trt — L relation only 


62 


56.3 


W a 


+ 2.0w 




0.6 ±0.3 


wo 


= ± 2 


60.3 


56.6 
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Fig. 1. — Lag - luminosity relation. Bursts with long lags are low luminosity, while bursts 
with short lags are high luminosity. The lag times for 39 GRBs have been corrected to 
the rest frame of the GRB and plotted versus the isotropic luminosity with the best fit 
power law (see equation 12) superposed. The one-sigma measurement uncertainties were 
used for the error bars. The lag-luminosity relation is thought to be caused by the speed 
of the radiative cooling of the shocked material (and hence the lag) being determined by 
the luminosity (Schaefer 2004). The observed slope in this Figure (-1.01) is close to the 
theoretically predicted slope of -1. The point of this Figure is that the lag time can be 
measured from 7-ray data alone and then a luminosity distance and a luminosity can be 
deduced for all GRBs. 
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Fig. 2. — Variability - luminosity relation. The variabilities for 51 GRBs have been corrected 
to the rest frame of the GRB and plotted versus the isotropic luminosity with the best fit 
power law (see equation 14) superposed. The one-sigma measurement uncertainties were 
used for the error bars. The variability-luminosity relation is thought to be caused by the 
fact that both the variability and luminosity are functions of the bulk Lorentz factor of the 
jet (Fjet)- The variability is a measure of the sharpness of the pulse structure, which is 
determined by the size of the visible region in the jet which is determined by Tj et . There is 
substantial scatter in this relation, so the variability carries less distance information than 
does the lag or E pea k- For low luminosity GRBs, the variability becomes small so that it 
is comparable or smaller than the ordinary Poisson noise in the light curve and hence the 
uncertainties in V become large. 
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Fig. 3. — - luminosity relation. As Tj et for a burst increases, the luminosity (L) 

increases as a power law of Tj et while the E pea k also increases proportionally to Tj et due to 
the blue shifting of the emitting region (Schaefer 2003b). This prediction that E pea k would 
be a luminosity indicator is confirmed by this figure. The E peak values for 64 GRBs have 
been corrected to the rest frame of the GRB and plotted versus the isotropic luminosity with 
the best fit power law (see equation 16) superposed. This relation is not the 'Amati relation' 
(the E peak -E 1:iso relation; Amati et al. 2002), which is based on different physics (similar to 
the E peak -E 1 relation). 



-77- 




Log[E peak (l+z)] 



Fig. 4. — E pea k - E 1 relation. The E pea k values for 27 GRBs have been corrected to the 
rest frame of the GRB and plotted versus the total burst energy in the 7-rays (i? 7 ) with the 
best fit power law (see equation 18) superposed. This relation is the tightest of the GRB 
luminosity indicators, and so the derived distance will carry higher weight than distances 
from, say, the V—L relation. To be included in this plot, the GRB afterglow must have 
an observed jet break in its light curve, and this means that only a fraction of GRBs with 
redshifts can yield a distance measure from this relation. This relation has been used by 
various groups (e.g., Ghirlanda, Ghisellini, & Lazzati 2004; Xu, Dai, & Liang 2005) but 
never with more than 18 GRBs included (Nava et al. 2006). This paper includes 9 more 
events (a 50% increase in the sample) mainly due to new bursts from Swift. 
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Fig. 5. — Minimum rise time - luminosity relation. This luminosity relation was originally 
predicted because the luminosity depends as a power law of the bulk Lorentz factor of the 
jet (Tjei) while the minimum rise time in the GRB light curve also depends on Tj et through 
the size of the visible region in the shocked jet (Schaefer 2002). This Figure provides a test 
and confirmation of the prediction that trt is a correlated with luminosity as a power law. 
The rise times for 62 GRBs have been corrected to the rest frame of the GRB and plotted 
versus the isotropic luminosity with the best fit power law (see equation 20) superposed. 
The one-sigma measurement uncertainties were used for the error bars. 
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Fig. 6. — N pea k - luminosity relation. With pulse rise times and durations depending on the 
opening angle of the visible region (and hence on Tj et ) , low luminosity bursts (with low Tj et ) 
will have very broad pulses which inevitably blur together no matter how many pulses are in 
the light curve. But a high luminosity event (with high Tj et ) will have short pulse durations 
so that pulses are unlikely to be overlapping. Thus, a theoretical prediction (Schaefer 2002) 
is that a lower limit on the burst luminosity can be derived from a count of the number of 
peaks in the light curve. This Figure shows a test and confirmation of this prediction, as 
all 64 GRBs lie above the theoretical limit (the curved line). To a close approximation for 
two-or-more peaks, the limit on luminosity is a power law (see equation 22) as displayed 
by a narrow line segment in the Figure. This relation only provides a lower limit on the 
luminosity, and thus is not of any use for constructing a GRB Hubble diagram. But this 
relation is useful for quickly putting limits on redshifts as an alert to observers, for example 
a faint Swift burst with many peaks is rapidly known to have high redshift. 
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Fig. 7.— GRB Hubble diagram for the concordance cosmology. For the concordance 
cosmology (for a flat universe with Q M = 0.27 and w = —1), this HD contains 69 GRBs 
out to redshift of >6. On this plot, look to see the redshift regime covered by supernova, 
with almost all being with z < 1 and only 10 with 1 < z < 1.7. The curve is the model as 
given by equation 10. We see that the observations follow smoothly along the model, indeed 
with chi-square of 72.3 and a reduced chi-square of 1.05. This is to say that the GRB HD is 
consistent with an unchanging Dark Energy. 
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Fig. 8. — GRB Hubble diagram for Riess' best fit cosmology. Riess et al. (2004) fitted 
a large sample of supernovae including those discovered by the Hubble Space Telescope 
(the so-called 'gold sample') and found that their best fit cosmology had the Dark Energy 
changing with time as described by wq = —1.31 and w' = 1.48. With this, the equation of 
state parameter for the Dark Energy (w) varies with redshift asw= —1.31 + lA8z. (This 
particular expansion has problems with divergences at z > 1 and will only be used in this 
paper for comparison with the Riess cosmology. Instead, the expansion of Linder (2003) will 
be adopted elsewhere.) Note that the observed points shift around by a small amount (with 
an RMS scatter of 0.10 mag about the average difference) as the luminosity relations have 
to be calibrated for every cosmology. We see that the bursts with z > 3 are systematically 
high, with only 4 out of 18 being below the model curve. The chi-square for Riess' cosmology 
is 80.3, which is substantially higher than the chi-square for the concordance cosmology. 
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Fig. 9. — Derived versus known redshifts. How good are the derived redshifts? For each 
burst, the luminosity indicators can be combined to give a distance modulus and then con- 
verted to a redshift for some adopted cosmology. This figure plots the derived versus the 
known redshifts for all 69 GRBs for the concordance cosmology. If the derived redshifts were 
perfect, then they would all be precisely along the diagonal line. The observed scatter about 
the diagonal has a reduced chi-square of 1.03, which indicates that the derived error bars are 
correct. The scatter about the diagonal does not vary with redshift, which means that the 
accuracy in the derived values are roughly some constant fraction. The RMS scatter about 
the diagonal shows that the average accuracy is 26%. 
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Fig. 10. — VLm for a flat universe with w — — 1. The best value and its asymmetric one-sigma 
uncertainty is 0.39to;Qg. 
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Fig. 11.— GRB HD results for evolving Dark Energy with a loose constraint on Q M . 
This plot shows the one-sigma and two-sigma confidence levels for a flat cosmology where 
Qm was constrained to be between and 1. The point of best fit implies a relatively 
shallow slope for z < 0.5 and this can already be ruled out with supernova data. Without 
marginalization, the best fit point has a chi-square of 70.7 for 66 degrees of freedom while 
the concordance cosmology has a chi-square of 72.3 for 69 degrees of freedom. This indicates 
that the concordance model is easily consistent with the GRB HD. 
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Fig. 12. — GRB HD results for evolving Dark Energy with a tight constraint on Q M - This 
plot shows the one-sigma and two-sigma confidence levels for a flat cosmology where M was 
constrained to be 0.27 ± 0.04. The difference between Figures 10 and 11 arises only from a 
differing constraint on Qm- The same analysis and arguments as in Figure 10 show that a 
constant Cosmological Constant is acceptable. 
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Fig. 13. — GRB HD with only the r /a9 — L, E peak — L, and E peak — relations. Do the 
V — L and t rt — L relations add more noise than signal? This 66 GRB HD was constructed 
without these two relations and for the concordance cosmology. This HD shows that the 
data still follow a smooth curve with no trends. The scatter about the model curve is larger 
than that in Figure 7, which is expected due to the V — L and Trt — L relations carrying 
32% of the statistical weight. This demonstrates that the combination of all five luminosity 
relations is greatly better than using any one or two luminosity relations. 
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Fig. 14. — GRB HD with no luminosity indicators. This small version of Figure 7 displays 
the HD for the case where the luminosity relations are taken to be flat (i.e., with their 
returned luminosities being a constant). This is a blatantly wrong assumption. The point 
of this figure is that this HD is horrible with huge scatter, and the comparison with Figure 
7 shows that the luminosity relations are producing a vast improvement. 
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Fig. 15. — HD for a representative range of cosmological models. The concordance cos- 
mology (with Qm = 0.27 in a flat universe with w = —1) is now the standard and default 
based on the supernova measures of the HD up to z < 1 with a handful up to 2 < 1.7. 
But many dozens of reasonable models for the equation of state for the Dark Energy as 
well as for alternatives to General Relativity have been proposed (cf. Szydlowski, Kurek, 
& Krawiec 2006). Three represntative alternatives are the Weyl Gravity (Mannheim 2006) 
with q = —0.2, a Chaplygin gas (Kamenschik, Moschella, & Pasquier 2001) with A=0.5, 
and the 'Riess cosmology' of w = — 1.31 + 1.482 as based on the best fit to the 'gold sample' 
of supernovae (Riess et al. 2004). Many of these models show miniscule differences up to 
z ~ 1.5 yet large differences for z > 3. This figure has three points: First, many reasonable 
alternatives to the concordance cosmology exist. Second, distinguishing between these cos- 
mologies at the relatively low redshifts available to supernovae will be hard due to the small 
differences that could get hidden by systematic errors (cf. Dommguez, Hoflich, & Straniero 
2001; 2003). Third, distinguishing between these cosmologies at redshifts > 3 is easy due to 
the large differences uniquely measured with GRBs. 



